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■ This is the second of a series of papers extending the 1+3 Covariant and 

Gauge-Invariant treatment of kinetic theory to an examination of Cosmic Mi- 
crowave Background temperature anisotropies arising from inhomogeneities 
in the early universe. The first paper (Paper I) dealt with algebraic issues, 
^ representing anisotropies in a covariant and gauge invariant way by means of 

^ , Projected Symmetric and Trace-Free tensors. 

Here we derive the mode form of the Integrated Boltzmann Equations, first, 
giving a covariant version of the standard derivation using the mode recur- 
sion relations, second, demonstrating the link to the the Multipolo Divergence 
Equations and finally various analytic ways of solving the resulting equations 
are discussed. 

A general integral form of solution is obtained for the equations with Thom- 
son scattering. The covariant Friedmann-Lemaitre multipole form of the trans- 
om ' port equations are found near tight-coupling using the covariant and gauge- 
invariant generalization of the Peebles and Yu expansion in Thompson scatter- 
ing time. The dispersion relations and damping scale are then obtained from 
the covariant approach. The equations are integrated to give the covariant 
^ ■ and gauge-invariant equivalent of the canonical scalar sourced anisotropics in 
^ ' the K = (fiat background) case. 

' We carry out a simple treatment of the matter dominated free-streaming 

^ ^ projection, slow decoupling, and tight-coupling cases in covariant and gauge- 

invariant theory, with the aim both giving a unified transparent derivation of 
this range of results and clarifying the formal connection between the usual ap- 
1 J proaches (for example Hu & Sugiyama) and the covariant and gauge-invariant 

' like treatments for scalar perturbations (for example of Challinor & Lasenby). 

1. INTRODUCTION 

The present series of papers (Paper I [18]) establishes the l-f 3 covariant kinetic theory for- 
mahsm of EUis, Treciokas and Matravers [29, 30] in a form that makes possible an investigation 
of the Cosmic Microwave Background anisotropies in the non-local context of emission of radi- 
ation at the surface of last scattering in the early universe, and its reception here and now (the 
'Sachs- Wolfe effect' and its later extensions). These papers aims to provide the link between 
the exact (non-linear) theory [59] and the linearised threading formalism (this paper and to 
a lesser degree [10]), to the linearised foliation formalism, based on Bardeens' approach [3] to 
cosmological perturbations [39] . 

The Ellis, Treciokas and Matravers papers introduced a covariant approach to kinetic theory 
based on a H-3 covariant representation of Cosmic Microwave Background anisotropies in 
terms of Projected Symmetric Trace-Free (PSTF) tensors orthogonal to a preferred time-like 
vector field m" [76, 67]. The benefits of the approach have been briefly summarized in Paper 
I [18] (see also Challinor & Lasenby [10, 11]) - we seek clarity by providing a direct formal 
derivation of the standard results as well as providing the background necessary to include 
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the non-pcrturbativc corrections discussed by Maartcns, Gcbbic and Ellis [59]. In essence, we 
provide (a) clear definitions of the variables used, (b) 1+3 Covariant and Gauge-Invariant (CGI) 
variables and equations, (c) a sound basis from which to proceed to non-linear calculations (as 
introduced in [59]), and (d) the possibility of using any desired coordinate and tetrad system 
for evaluating the variables and solving the equations in specific cases (because the general 
equations and variables used are covariant). 

This approach has been used in a previous series of papers [57, 73, 58] to look at the local 
generation of anisotropies in freely-propagating radiation caused by anisotropies and inhomo- 
geneities in any universe model, providing a proof that near-spatial homogeneity in a region 
U follows from radiation ncar-isotropy in that region (an 'Almost-Ehlers, Geren and Sachs' 
theorem, generalizing the important paper by Ehlers, Geren, and Sachs [23]). By contrast, 
the present series looks at non-local (integrated) anisotropy effects in the context of the stan- 
dard model of cosmology - given the observational justification provided by the almost-Ehlers, 
Geren and Sachs theorem [72]. 

Paper I [18] of the present series of papers [18]) dealt with the CGI irreducible representation 
of cosmic background radiation anisotropics by PSTF tensors, and their relation to observable 
quantities (specifically, the angular correlation functions) , first within a general framework and 
then and specialized to almost-Priedmann-Lemaitre universe models ^ [57], as well as dealing 
with multipole and mode expansions and the relation to the usual formalisms in the literature 
[82, 35, 39, 40, 20]. 

In this paper (Paper II), we use the CGI approach to study cosmic background radiation 

anisotropies in almost- Fricdmann-Lcmaitre models in an analytic manner, by time-like integra- 
tion of the almost-Friedmann-Lemaitre differential relations. Our emphasis is on the canonical 
linearised model for cosmic background radiation anisotropies [68, 81, 1, 75, 35, 39, 40, 79], 
systematically developing the CGI approach and providing a comprehensive and transparent 
analytic link to the alternative analytic gauge- invariant treatments based on the Bardeen gauge- 
invariant variables. We develop these results both in mode and multipole form, emphasizing 
the different physical processes and assumptions and demonstrating how these are dealt with in 
the CGI context. This requires the covariant mode form of the integrated Boltzmann equations, 
based on the recursion relations for almost-Friedmann-Lemaitre mode functions [18], enabling 
a direct mirroring of standard treatments based on Wilson's seminal paper [81, 35, 39, 40, 79], 
except carried out in a CGI fashion, thus forming a sound basis for extension to non-linear 
effects. 

Indeed one of the advantages of the CGI approach in the context of the generic multipole 
divergence equations is the ability to include non-linear corrections to the almost-Friedmann- 
Lemaitre treatment. Towards this end, the relationship between the covariant mode formula- 
tions and the almost-Friedmann-Lemaitre covariant multipole treatment arc given with this in 
mind, based on the relations in Paper I [18]. A key point here is that there is nothing new 
about the linear formulation itself, however recovering the standard analytic linear results from 
the the exact theory [59] in a straightforward way, is new. 

Moreover, the results presented here provide the foundation for a non-linear extension of this 
approach which is outlined in a paper by Maartens, Gebbie and Ellis [59]. 

We emphasize that in our treatment, (t^^t ) (the multipole form of the angular correlation 
function) is given for small temperature anisotropies irrespective of the form of the geometry 
[59], making it the natural representation for the inclusion of non-linear dynamic effects, while 
the analysis for [t^I^ (the mode coefficient mean square) is specifically for almost-Friedmann- 
Lemaitre models [18]. The non-linear extension of the almost-Friedmann-Lemaitre treatment 
given here will be based on the multipole-to-mode relations, leaving the use of mode expansion 
to the latest possible stage. 



^In order to be clear on the use of these names, Robertson- Walker refers to the Robertson- Walker geom- 
etry whatever the dynamics, while Priedmann-Lemaitre refers to such a geometry which additionally obeys 
the Priedmann-Lemaitre dynamics implied by imposing the Einstein Field Equations. An almost-Friedmann- 
Lemaitre universe is a universe model governed by the Einstein field equations, whose difference from a 
Priedmann-Lemaitre universe is at most 0(e) in terms of a small parameter e [25]. 
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Our focus is on the era following spectral decoupling (near 500 eV). A complex series of 
interactions take place at the various epochs of the expansion of the universe. The kinetic 
equations developed in Paper I and [59] can represent almost any such interactions provided we 
use appropriate interaction ('collision') terms; the issue is how to obtain simplified models that 
are reasonably accurate in the various epochs. We will consider only two kinds of representation 
here: namely 

1. Thomson scattering, valid at late times when particle and photon numbers are conserved 
and the matter is non-relativistic (during decoupling, an alternative approach is to use a 
visibility function). 

2. An effective two-fluid approach, obtained by truncation of the hierarchy and valid at 
earlier times when strong interactions take place establishing equilibrium or close to equilibrium 

conditions between the components, i.e. when the interaction time-scale is much less than the 
expansion time-scale; an alternative description is to use a single imperfect fluid [60]. 

Both descriptions are valid when the matter is relativistic. The detailed form of interactions 
does not need to be represented in this case, because the state of the matter depends only on 
its equilibrium condition, characterized by its equation of state. 

At some times either form is valid and they can then be related to each other. We do 
not attempt here to give a description of earlier non-equilibrium eras when processes such as 
pair production, nucleosynthesis, etc, occur, nor do we consider issues such as inflation and 
reheating after inflation, and the differences between the inflation sourced perturbations as 
opposed to those based on other phase transitions. Thus our models will be valid after the end 
of any period of inflation that may have occurred and after strong non-equilibrium processes 
have ceased. During this era the processes which result from initial fluctuations left over from 
earlier non-equilibrium epochs, determine the final cosmic background radiation anisotropy. 

Specifically, we deal with four eras of interest. Going backwards in time from the present, 
they are, firstly, free streaming from last scattering to here and now, in the matter dominated 
almost-Friedmann-Lemaitre context; secondly, slow decoupling during a matter dominated 
era, which is when the cosmic background radiation spectrum freezes out; thirdly, the late 
tight-coupling era after matter-radiation equality, during which structure formation begins; 
and fourthly, the main tight coupling era after any inflationary epoch and before matter- 
radiation equality, during which acoustic modes occur in the tightly-coupled fluid, the initial 
matter perturbations having been seeded by earlier conditions (for example, inflation). We then 
show how to put these CGI results together to determine the major features of the expected 
anisotropy spectrum. We develop sound models of the dominant effects in each of the eras we 
consider, but there will always be a ncxxl for reflnement of these models by taking into account 
further effects (in particular polarisation and the effect of neutrinos). 

Although the effect of the neutrinos is crucial, and can been subsumed into the gravitational 
variables. We do not provide nor discuss further the additional hierarchy of neutrino moment 
equations that would then need to be included. The exact massless neutrino evolution equations 
are given in a previous paper [59] and it is easy to show that the resulting linearised equations 
are essentially the same as those for massless radiation without coUisionals, hence in this paper 
we focus on the photon equations. 

In more detail: 

• Free- Streaming : We find the CGI integral solutions to the almost-Priedmann-Lemaitre 
multipole divergence equations with no collision term, and use them to project the initial data 
from decoupling to the current sky. We explicitly do this for instant decoupling. Neither 
the Vishniac, Rees-Sciama, thermal Sunyaev-Zel'dovich nor lensing effects are considered here 
— we will focus on the CGI model of the dominant proc(«sc;s, and these hirther effects will 
introduce detailed modifications. However, a comprehensive understanding of such secondary 
higher order effects relies on a derivation of the anisotropy effects given here. 

• Slow- Decoupling : Here we consider modification of the previous results when slow decou- 
pling of the interactions due to Thomson scattering is taken into account. We consider the 
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damped integral solution for slow recombination, and as an alternative description modify the 
integral solutions appropriately with a visibility function carrying the functional dependence 
of the varying electron fraction, in a matter dominated context. Effectively, recombination is 
complete before the radiation decouples. This means that the surface of last scattering is found 
a little after the end of recombination ^. It is during this era that photon diffusion damping 
scale effects become important - the damping scale is affected by the duration of this era. 
This will be investigated in the context of almost-Priedmann-Lemaitre universes after matter 
radiation equality. 

• Tight-coupling: This is the key to the entire treatment. We give the CGI version of 

the tight-coupling approximation of Hu et al [39, 40, 44]. In the almost-Friedmann-Lemaitre 
treatment, the slow decoupling and free-streaming era's will only 'damp' and 'project' the 
spectrum formed at the end of tight-coupling onto the current sky ^. 

We consider two different tight-coupling regimes. 

— The late tight coupling era is separated, conceptually, from the early tight coupling 
era by matter-radiation equality, after which time the matter perturbations have effectively 
decoupled and (CDM-based) structure formation begins. In this era, strong interactions have 
ceased and a Thomson scattering description can be used. We first carry out the near tight- 
coupling treatment of this era based on Peebles & Yu [66, 47], and then reduce these equations 
to the covariant tight-coupling equations equivalent to those of Hu et al [39, 40, 44]. This 
provides the basis for understanding the acoustic signatures in the temperature anisotropics 
within the CGI approach. 

— The early tight coupling era occurs between the matter and radiation eras (after strong 
interactions have ceased), when a Thomson scattering description will also be sufHcient. This 
era is characterized by acoustic oscillations in a tightly coupled fluid; for calculation convenience 
this can be represented as a single dissipative fluid [60], or for a slightly more sophisticated 
treatment by tightly coupled two- fluid models [14, 13]. We give a CGI derivation of the 
harmonic oscillator equation providing the primary source terms in the standard model of 
Doppler peak formation by acoustic oscillations. 

We put these results together in sections 7 and 8, where the equations are integrated to give 
the CGI treatment in the K = (flat background) case in terms of an integral solution. 

The primary sources of the temperature anisotropics (the acoustic and Doppler contributions 
near last scattering resulting in 'Doppler Peaks' today) arc demonstrated ^. This recovers the 
Sachs- Wolfe family of effects for flat background Robertson- Walker geometries, but derived 
from a CGI kinetic theory viewpoint as opposed to the photon-propagation description used 
in the original Sachs- Wolfe paper. 

The form of the angular correlation function is determined for the primary effects (although 
not given explicitly in terms of the matter power spectrum). The normalization to standard 
CDM (Adiabatic CDM) is presented in terms of CGI variables. This demonstrates the basic 
effects in the CGI formalism, and links our approach to the standard literature, see for example 
[79] and references therein, where further details of this era are given. 

It is important to note that the integrations considered here are carried out along time-like 
curves, even though the cosmic background radiation reaches us along null curves. These are 
alternative approaches that are equivalent in the context of linearisations about Priedmann- 
Lemaitre models; differences will however occur if we include non-linear corrections. Briefly, 
the key point about cosmic background radiation integrations is that there arc two ways in 

note on nomenclature: By decoupling we have in mind the situation when the interaction rate per particle 
becomes less that the expansion rate. By last scattering surface we mean the surface upon which the diffusion 
scale is equal to the horizon scale, after which it is larger than that scale and the free-streaming approximation 
is sufficient. The photons will decouple from the thermal plasma near 0.2 eV, and from the matter after 
recombination has effectively ended, near 0.3 eV. Free-streaming is considered to be a good approximation from 
about 0.26 eV. 

^This useful feature of the almost-Priedmann-Lemaitre models is due to the homogeneity and isotropy con- 
ditions in the background Priedmann-Lemaitre universe and is not generic [59]. 
■^following the analytic treatment of Hu & Sugiyama [39, 40] 
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which to proceed: Firstly a null-cone integration, following the radiation from last scattering to 
the present day^, and secondly a time-like integration along the matter flow lines (as here). In 
the latter case one is (at least implicitly) thinking of a small comoving box containing matter 
and radiation [36] which is similar to all other small boxes at the same time, and where one 
has assumed that the radiation leaving is exactly balanced by the radiation entering (from 
neighbouring boxes), whether in tight-coupling (when it is a local assumption) or in the free- 
streaming era (when it is a non-local assumption) . In effect one integrates behaviour in such a 
box in a small domain about our own world-line from tight-coupling through decoupling to the 
present day; to do this, one does not need to know about the behaviour of null-geodesics (the 
integration is along time- like geodesies). Before decoupling the matter and radiation evolve 
as a unit while after they need to be integrated separately (giving the corresponding transfer 
functions in each case). Then this is related to observations by, first, conceptually shifting 
copies of the small box at the time of last scattering from our world line to all points where the 
past null-cone intersects the surface of last scattering at that time; this can be done because 
spatial homogeneity says that these boxes are essentially the same (a Copernican assumption 
is used here, justified by the almost-Ehlers-Geren-Sachs theorem [57, 58]); second, by then 
relating distances on the last scattering surface to observed angles by using the area distance 
relation, relating physical distances at last scattering to angular size in the sky. 

The next paper in the series. Paper III [15], deals with the explicit relationship between 
null-cone and time-like integrations. Further papers will look at non-linear extensions of the 
results given here. 



2. LINEARISED COVARIANT MODE EQUATIONS 

To study details of cosmic background radiation anisotropy generation we need both a spa- 
tial Fourier decomposition, defining wavelengths of perturbations, together with the angular 
harmonic decomposition relating anisotropics to observed angles in the sky. The CGI versions 
of both decompositions were given in Paper I [18], giving the relationship between the mode 
and multipole variables. 

The dynamic relations obeyed by these quantities, determining the cosmic background ra- 
diation spatial and angular structure, can be obtained from the Boltzmann equation in two 
ways: via multipole divergence equations or via the integrated Boltzmann equations. In each 
case the general equation needs to be mode-analyzed to obtain the spatial structure. 

In the first case, the almost-Friedmann-Lemaitre multipole divergence equations are obtained 
by systematically linearising the non-linear multipole divergence equations for small tempera- 
ture anisotropics given in Maartens, Gebbie and Ellis [59] . The mode form of these equations 
(16-18) [59] can then be obtained by mode analysis, see (19-20) below. 

By contrast, the more common procedure is to directly construct the mode form of the 
integrated Boltzmann equations from the linearised integrated Boltzmann equations by sub- 
stituting (B.3) into (B.l) and integrating over the energy shell with respect to E^dE (For a 
more detailed relativistic kinetic theory description of these equations see [59] and [11]): 



f 

Jo 



r 1 8fi r°° 

E^dE [e{u'' + e-) Va/ - ( 3 e + A„e° + aabe''e'')E^ « E^dEC[f] . (1) 



Upon using the CGI definition of directional bolometric brightness [59] : 



T{x) [1 + T(x,e)] = 



47r 
r 



j E^f{x\E,e'')dE 



1/4 



(2) 



the covariant equivalent of the standard formulae given in the Wilson-Silk approach [81] can 
be found, where as usual the partial energy derivative is removed by a integration by parts and 



^which can be parametrized either by a null cone parameter, a projected spatial coordinate, or a projected 
time coordinate 
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application of the regularity conditions. This is the approach we develop now, showing later 
the relation to the approach based on the multipole divergence equations. 

2.1. The Mode Equations 

The optical depth k is given in terms of the Thompson scattering cross-section ar, the 
number density of electrons rig, and the fraction of electrons ionized Xg- 

K{to,t) = aTne{t')xe{t')dt' = J kdt' , (3) 

where r] is the conformal time {dt = adrj) in the Friedmann-Lemaitre background. Starting 
with the almost-Friedmann-Lemaitre integrated Boltzmann equations (1) for the temperature 
anisotropics (2) with the coUisional term for isotropic scattering in terms of the optical depth, 
and the expansion replaced by substituting from the radiation energy conservation equation 
(G.11,15), we obtain the linearised integrated Boltzmann equations: 

-r « e»D„r - iD„r° + (D„ InT + A„)e« + aabC^e^ - K(e«t;f - r) . (4) 

We take a mode expansion (see Paper I [18]) for r (the temperature anisotropy), Aa (the ac- 
celeration), DQ(lnT) (the spatial-temperature perturbation), (Jab (the shear), and the gradient 
of the radiation dipole, D"ra, based on solutions Q{x) of the scalar Helmholtz equation: 

D«D„Q = -%Q (5) 

in the (background) space sections, where the Q's are covariantly constant scalar functions (i.e. 
to linear order Q w 0) corresponding to a wavenumber k. The physical wavenumber is defined 
by kphys{t) = k/a{t). These functions define tensors QAi{k'^,x'^) that are PSTF (in the case 
of scalar perturbations they are given by the PSTF covariant derivatives of the eigenfunctions 
Q) and so allow us to define the mode functions [18]: 

QAe = i--] ^{A,)Q and Ge[Q] = 0^'QAe , (6) 

where O^* = e^'^*'* is the trace-free part of We can expand any given function f{x,e) in 
terms of these functions, thus for scalar perturbations (see (G.20 -G.23), the mode coefficients 
of the temperature anisotropy are constructed as follows: 

oo 

T(a;,e)=^^Tf(i,fc)G4g]. (7) 
e=i k 

Note that tq is identically zero, because (2) defines the temperature T gauge-invariantly as 
the all-sky average in the real (lumpy) universe (it is not defined in terms of a background 
model). We can then write: 

DV„ = ^-Ti(fc,t)Q, Ba{lnT) + Aa=Y^[liST{k,t)+A{k,t)]Qa, (8) 
k k 

and 

<yab = ^CF{k,t)Qab- (9) 

k 

The mode coefficient of the velocity of the baryons vb relative to the reference frame is 
given by 
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The equations arc linearised at 0(?j^), 0(ei;) and O(e^) [59], so we can use, for example, 
u% u"" + v% to give the baryon relative velocity. 

The equations here are gauge-invariant (given relative to a unique physically-based choice 
of the 4- velocity vector m") and valid for any choice of mode functions, but the detailed result 
of their translation back into the spacetime representation t{x, e) via (7) will depend on the 
harmonic functions Q{x^) chosen^. 

Wc substitute these expansions into the linearised integrated Boltzmann equations and then 
use the recursion relation [81, 35, 39, 18] for mode functions G£[Q]{e°' , x^) with wave-number 
k : 



e^D^lGelQ]] = - 



(2f -M)(2£- 1) 



1) G^_i[g]-G^+i[0] 



(10) 



where K is the curvature constant of the background space sections. 

With the mode decompositions of each term in (4) for each wave number^, on using the 
recursion relation (10) and separating out the different harmonic components, the almost- 
Priedmann-Lemaitre mode equations are found from (3)^: 



-re 



-T2 



-Tl 



(2^ + 3) (2^+1) 



+ kTe ,i>3, 



All 

35 



8K 



T3 - Tl 



KT2, 



15V^"¥'^^ 



[^ST + A]-k{vB-n) 



(11) 
(12) 

(13) 



The above equations demonstrate the up and down cascading effect whereby lower order terms 
generate anisotropics in the higher order terms, and vice versa, in a wavelength-dependent way; 
curvature affects the down-cascade but not the up one. These equations can be compared to 
the equations of Hu & Sugiyama, in particular (eqn. 6, p. 2601) [40]. They are identical if we 
use the Newtonian frame (discussed in the following sections), and so have the same physical 
content, but are more general since they are valid with respect to a general frame m". 

2.2. Prom Multipole Equations to Mode Equations 

The relationship between the angular harmonic and mode expansions are given in Paper I 
[18]. We start by writing the CGI harmonic coefficients ta^ in terms of the mode functions 
(6): 



TAe = Y,Te{k,t)QAe « J2^eit,k){--)-^D^A,)Q 



(14) 



Then the angular harmonic expansion for / becomes the mode expansion (7). On taking 
the multipole integrals of / as in Paper I [18], they too are then mode-expanded by (14); so 
the linearised divergence relations for these multipoles given in [59] become mode equations, 
equivalent to the almost-Friedmann-Lemaitre mode equations (11-13) derived above. 
In detail: The almost-Priedmann-Lemaitre multipole divergence equations are [59] : 



(15) 




®In effect there are two major choices, namely plane wave solutions and spherical solutions; the former occur 
naturally in galaxy formation studies and the latter in observational analysis, so the relation between the two 
(see Paper I [18]) is a central feature of analyzing null-cone observations. 

'^Tlicrc should be a summation over wave numbers in the following equations. However we follow the estab- 
lished custom in omitting this summation and any explicit reference to the assumed wave number k. 

®Note that these equations are valid for any choice of Q, including both spherical and plane wave harmonic 
functions. 
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i-fa) ^ T>alnT + Aa + p^Tac - CTTUeiv^ - Ta), (16) 
{-fab) - CTab + D(oT"6) + f D'^Tobc + (^TTleTab, (17) 

i-TAe) ^ D(„^T^^_^) + ^^i-^DVA.c + O-TneTAf ■ (18) 



Now the following identities are used (dropping the k-summation): 

' O^^Qa, , (19) 



0^'B^,^TA,_,) « T,_i ( ~ ] O^^QAe , (20) 



where the first relation arises from the use of the identity 



(£ + 1) ( k 



(2£+ 1) V a 



U2 



-,2 



1- 1^(^ + 2) 



These are substituted directly into the multipole equations after taking a mode expansion of 
those equations and then dropping the k-space summation, loading again to the equations 
(13). The point to note is that while one does not explicitly need the multipole equations 
in order to find the almost-Friedmann-Lemaitre mode equations (which can be derived from 
the linearised integrated Boltzmann equations as shown above), in order to examine non- 
linear effects one can obtain the necessary equations by proceeding as here from the non-linear 
multipole divergence equations, to obtain higher approximations of the mode equations and 
the mode- mode couplings. 

2.3. The Einstein Equations 

The key quantities which link the radiation evolution through to the matter in the spacetime 

geometry arc the shear aab = U{a;b): the acceleration Aa = Ua-.bU^ and the expansion 0. These 
couple the multipole divergence equations to the Einstein field equations (which are given in 
Appendix G, see (G.3)-(G.18)). 

The shear and acceleration, arc related to the electric part of the Weyl tensor -E(ab), the 
anisotropic pressure iri^ab) and matter spatial gradients (see (G.3) -(G.18)) while the CGI spatial 
gradient of the expansion is linked to the divergence of the shear, heat flux vector Qa and the 
vorticity vector Ua- 

-^{p + p)aab « {Eab + hi^ab)+3H{Eab + ^Trab)-HTrab-{h^{aqb)}, (21) 

{p+p)Aa « -Dap-'D\ab-{qa+4.Hqa}, (22) 

iD„e « liB'aab) -{ya + curl a;,} . (23) 

These equations are valid for general almost-Friedmann-Lemaitre perturbations. In the re- 
stricted case of scalar perturbations, we set the vorticity to zero^'' and non-zero quantities can 
be written in terms of potentials [71]. In particular 

Eab ~ D(„Db)$B = iD(,Db)($A - $ff), (24) 
TTab « D(„Df,)$^ = -D(„Db)($ff + 3>a) , (25) 



^This has also been derived by Challinor and Lasenby [11] and is found from the PSTF recursion relations 
and the generalized Helmholtz equation (which is in turn found from the constant curvature Ricci identity) 
[18]. 

'^"We can obtain scalar equations even when the vorticity is not zero, by talking the total divergence of these 
equations; we will not pursue that case here. 
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where the potentials and >S/f arc; analogous to the GI potentials used by Bardeen [8]. The 
following useful combinations can be found : 

Eab- ^nab^'D{aDb)^A, and Kb + ^TTab « -D(„D(,)$h- (26) 

Using the Einstein field equations, the total flux, qa, can also be expressed covariantly in terms 
of these potentials: 

Hqa « B^D^^Dt,)^H - lT)aP, (27) 
HB^aQb) « D^.D,) [|(D2$H + (p-3/f2)$^)-ip] . (28) 

This then allows us to write the shear and acceleration in terms of the scalar potentials and 
perturbation variables: 

UP+PWb « (D^.D,)*^) +3HD(„D5)$H-i/D(„Db)($H + $A) + {iD(„g6)} , (29) 

(p + p)A, « -D„p - D^D<„Db) i^H + ^a) - {qa + ^^Hqa} . (30) 



2.4. Frame Transformations and Gauge Fixing 

There is freedom associated with the choice of rcJercmce velocity m", which we call a frame 
choice. This is to be distinguished from the choice of coordinates in the realistic universe model, 
which can be done independently of the choice of It is equivalent to the choice of time- like 
world-lines mapped into each other by the perturbation gauge chosen (i.e. the mapping between 
the background model and the realistic lumpy universe model, see for example Bruni and Ellis 
[25]), but is independent of the choice of time surfaces in that mapping. Given a particular 
covariantly defined choice for this velocity, the frame choice is physically specified and the 
equations are covariantly determined and gauge invariant under the remaining gauge freedom 
In simple situations this choice will be unique, however in more complex situations several 
choices of this velocity are possible, each leading to a somewhat different CGI description. 

When we restrict ourselves to a particular frame in order to simplify calculations, we can 
straightforwardly make the appropriate simplifications in the general equations to see what the 
implications are (for example setting = for the energy frame, the quantities in the braces 
in equations (29,30) above vanish). However it is also useful to explicitly transform from one 
frame to another and examine the resulting effect on dynamic and kinematic quantities. 

Under a frame transformation ps u°- + v"", \v°'\ <C 1 (restricting our analysis to non- 
relativistic relative velocities ^^), the following relations [59] hold: 





a aab + D(o^^6)) 


(31) 




a Aa + Va+ HVa, 


(32) 


e ? 




(33) 




- Qa- {P + P)Va, 


(34) 




a - IcmlVa . 


(35) 



The quantities p, p, TTab, Eab and Hab, remain unchanged to linear order in almost Friedmann- 
Lemaitre universes (e.g. Wab ~ T^ab and Eab ~ Eab), and the temperature anisotropies (r^^,) for 

-"-^ Gauge fixing requires in addition a sjjccification of the eorrespondcnce of time surfaces in the realistic and 
background models (eff'ectively specified by determining the choice of surfaces of constant time in the realistic 
universe model) and of points in initial space-like surfaces. 

^^It is important to recall that gauge invariance is only guaranteed if the choice of velocity Ua coincides 
exactly with its value in the background spacetime. This is not difficult in practice, as appropriate physically 
defined frames Ua will necessajrily obey this condition because of the Robertson- Walker symmetries. 
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£ > 1 arc similarly invariant for the small velocity transformations. The baryon and radiation 
(dipole) relative velocities change according to: 

Va , (36) 
Va ■ (37) 

Also, the projection tensor hab changes if we boost to the frame Ua giving hat, hence any 
spatial gradients need to be modified and will be the totally projected derivative in that 
frame. The consequence is that the perturbation variables change accordingly: Thus for any 
species I, [59] 

Bain PI « Balnpi - 3iJt;a(/ + /)// • (38) 

For example I = R and I = B give the equations for radiation and baryons respectively, 
implying: 

BalnT ^BalnT -VaH , D„ In pm « D„ In pm - 3ff w„ . (39) 




These equations allow us to determine the required transformation to obtain desired properties 
of a particular choice m". The almost-Friedmann-Lemaitre multipole divergence equations (16- 
15) are valid in any frame; in particular, if a frame transformation is performed as above, they 
can be given in terms of the resulting variables in the new frame, Ua, with whatever restrictions 
result. 

While various choices of {t" are available in a multi-fluid description of the early universe 
[59], there are three particularly useful choices. 

1. The energy frame: = is preferred when dealing with two coupled particle species, 
as in the two fluid scenario [13]. This is useful as the Einstein field equations are simplified 
to a form which takes on a similar structure to the matter dominated equations, and even in 
the strong interaction case may be expected to be unaffected by collisions because of energy- 
momentum conservation (this choice is dealt with in more detail below in the context of scalar 
perturbations). 

2. The zero acceleration frame (or CDM frame): Ua = ~ Ua -\- v'^ . Prom the CDM 
velocity equation [11, 59] and (32) we then find: i-f + + Aa ~ Q ^ ~ [56, 11, 59]. 
This choice is particularly useful in multi-species situations, as this frame will be geodesic right 
through tight-coupling, slow decoupling and into the free-streaming era. 

3. The Newtonian frame: Ua = Ua in which the vorticity and shear of the reference frame 
vanishes: dab ~ Bi^a^b) = 0, when such a frame can be found. This frame is only consistent in 
restricted cases [31] , but is particularly useful in making comparisons with much of the analytic 
literature [39, 40, 35] and in making connections with the local physics in terms of Newtonian 
analogues in Eulerian coordinates. For example, the matter shear can be then thought of in 
terms of distortion due to the relative velocities (31): aab ~ ~^{aV^)- 

4. The constant expansion frame: Do© = 0. This choice is sometimes useful when discussing 
perturbations on small scales. 

These various choices will simplify the equations in significant ways, and enable us to recover 
many of the standard results. It should be noted however, that the covariant equations we 
have given above are general and do not require either gauge or coordinate restrictions to be 
meaningful, and the covariant quantities have in them a natural invariant geometric meaning. 
We will therefore retain the covariant form and do not restrict ourselves to any particular frame 
nor gauge choice for most of this paper, however we retain the freedom to make such a choice 
when useful. If and when we do pick a particular frame, this will be explicitly stated along 
with the reason for doing so. 

2.4..!. The Newtonian Frame Link to the Bardeen Variables 
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Here we demonstrate the direct link between our variables in the scalar case, and those used 
in the Newtonian gauge, in terms of the Bardeen variables. From (31) and (32) we find easily 
that for u°' = where D^^nh) = 0, Ha = Ua + (the consistency of this choice is discussed 
in [31]): 



f 




(40) 


Aa ^ 




(41) 


DalnT P 




(42) 






(43) 


e p 


a e + divt;^. 


(44) 


Qa - 


- qa + {p + p)v^. 


(45) 



The efi'ect of this frame transformation is to modify the £ = I and 2 multipole divergence 
equations (15-17): 

-fa « DalnT + Aa + lB^Tab-aTn,{v^ -Ta), (46) 
-fab « D^„f(,) + ^B^Tabc + (TTTleTab ■ (47) 

The £ > 2 equations (18) remain unchanged, however the field equations as well the perturba- 
tion equations need to modified, if necessary using the transformation relations (32 -37). For 
example, (15) becomes 

{DalnT) + HiBalnT + Aa) ~ -iD„e - iD„(D-re), (48) 

which can easily be checked to be invariant under the frame transformations Ua ^ Ua + Va- 
From the shear evolution equation (G.6) : 

D(„ib) W Eab - ^TTab, ia « Da^A , (49) 

the momentum constraint (G.7) and the propagation equation for the electric part of the Weyl 
tensor (G.12) one finds respectively that: 

iD(,Db)e « -^B^aQb), |D(,D„)e « D(„D;,)$ff - HB^a^.^^A ■ (50) 

This gives the scalar monopole equation for the temperature perturbation: 

D(„D,>(lnT) « -D^„Db>$H - iD(„D(,>(DVe) • (51) 

It is important to note here is that in the shear-free frame we can interpret the acceleration 
directly in terms of the $a potential, in other words, in terms of its Newtonian analog, while 
can be interpreted as a curvature perturbation. In terms of the potentials used by [39, 16] 
one can identify \E' = and $ = 

The above formulation is useful in linking the covariant work to the usual GI treatments. So 
for example we take the mode expansion of the potentials, one finds on dropping the k-index 
on the right, 

Ac « ($A|a + + HVs^a) = {V^^'^ + ^V^f - k^A)Y('\ (52) 
ac,0 « a(D<„D^)ys) = -akV^°^Y^°\ (53) 
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where the prime (') denotes the time-derivative with respect to the conformal time, Y arc 
the eigenfunctions of Y^"\a = —k'^Y and following Kodama and Sasaki [50, 8], the bar (i^,) 
denotes spatial derivatives with respect to surfaces of constant curvature in the background. 
Furthermore, one can identiiy Vs as a relative velocity. 

2.4-2. The Energy Frame 

In order to be clear on the consequences and subtleties involved in fixing the frame, here we 
give the source terms in the energy frame (u" = =^ 9a = 0) for scalar perturbations. The 
important point is that this is a physical frame, uniquely defined by the local physics. The 
equations (29) and (30) then take on the form : 

{p + p)^ab « 2(D(„D6)$H)' + 4iJD(„Db)$H-2i?D(„Db)$A + D(„t;b^, (54) 

{p + p)A, « -D„p-D*D<„D(,>($H + $A)+i'f + (55) 

Many CGI treatments use this frame [13], and have the advantage that the equations take on 
a form which is similar to those for the matter dominated case, but can still be used near to 
matter-radiation equality. 

2.4-3. Matter Domination 

During matter domination (we have in mind the CDM dominated case) there is a unique 
physically relevant frame defined by the matter 4- velocity, hence without restricting the almost- 
Priedmann-Lemaitre universe further there is natural frame, u", in which the variables will be 
gauge invariant and the above relations hold. 

In this frame, equations (29) and (30) become 

a^p^aab « -2(a^D(„Db)$ijy, and A„ « . (56) 

Here p ~ p„ is now the density of the matter content only. The key point is that to retain 
a consistent linearisation scheme as well as retaining gauge invariance, we now have three 
smallness parameters: v (non-rclativistic relative velocities), 77 (radiation-baryon ratio is at 
least 10~^), e (the universe is almost Friedmann-Lemaitre when e is at least 10~^) [57, 59, 72]. 
It follows that (the radiation energy density) is now 0{r]) and we can drop all terms at least 
0{rie), O(e^) and 0{ri^) such as, for example. p<7ab ~ or ^p„Tab — T^ab ~ 0. This is how the 
anisotropic pressure is eliminated to the order of the calculation in this scheme. 

The link to the matter distribution in the spacetime comes through the mode coefficients 
a{k,t) (of the shear), A{k,t) (of the acceleration), and 5T{k,t) (of the temperature perturba- 
tion) . A mode analysis leads to a particular solution of the linearised Einstein field equations 
due to scalar modes as in (G.20 - G.23): 

Eab « ^Qab = '^^H{t, k)Qab, (57) 

(Tab- -I {H^^o)- ^e{a<^H(kMQab, (58) 
D„lnp„ « lk{Hl^^)-^-^H{k,t) [e-iK]Q,. (59) 

For a matter dominated open model {K ^ 0) where oq = +1 we have Hq k, K/{^q — 1). 
If we add the adiabatic assumption (see Appendix G.4) we find 

D„lnT« +iD„ln/)M, (60) 

where we have used that w 3HQ^lQa~^ in the background. One can then put the mode 
coefficients, in the matter dominated scalar adiabatic almost-Friedmann-Lemaitre models, into 
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the form : 

ST{k,t)^^{H^n,)-\a^H) [lik^-SK)] , (61) 

A{k,t)f^Ofiiv^{k,t), (62) 

aik,t) « ~^iH^Qo)-\a^H) [k^] • (63) 

The first expression gives the direct effect of the gravitational potential on the cosmic back- 
ground radiation anisotropies (Sachs-Wolfe effect), and the third the effect of the time variation 
of the potential on these anisotropies. These are investigated in detail in later sections. The 
matter dominated Einstein field equations are at least 0{eri) and fix the form of the shear, 
the acceleration and the temperature perturbations Da In T as they enter the integrated Boltz- 
mann equations (which is how the geometry enters into these equations). The hierarchy itself is 
0(e) and although the radiation variables do not enter the almost- Friedmann-Lemaitre (mat- 
ter) Einstein field equations, they remain non-zero, and therefore there are still temperature 
anisotropies, ta^- This is an important but subtle point - matter domination implies the 
radiation moves as a test field over the geometry. 

2.5. Linearisations, Approximations and Scales 

In this section we discuss the various linearisation and approximations (already mentioned 
in the last section), that will be used in this paper. 

2.5.1. Almost- Friedmann-Lemaitre Linearisation 

Here we drop all terms that are at least C(e^). The implication of this is that one can 
only consider small velocity boosts, large ones would break the linearity about the Friedmann- 
Lemaitre background - hence we include v'^ = \v"'Va\ <C 1 as an almost-Friedmann-Lemaitre 
limit dropping the additional terms that are at least 0{ve,v'^). 

The important subtlety here is that ef is in fact the temperature anisotropy smallncss param- 
eter as related to the temperature moment mean squares ItaJ oc ei. The almost-Friedmann- 
Lemaitre limits on the geometry, e, (which define Uab, Aa and -Da© (for example) as 0(e) in 
appropriate dimensionless units [73, 72] are related to ee via the almost-Ehlers-Geren-Sachs 
theorem. In other words, limits on the temperature anisotropics, e^, put bounds on the size 
of the smallness parameter e, given that a weak Copernican principle holds. Furthermore, the 
limits on e in turn place consistency limits on the size of the v/c boosts that are applicable 
(here in units of c = 1). Thus at least almost-Friedmann-Lemaitre means keeping terms that 
are at most : 

Almost-Friedmann-Lemaitre « 0{e,v). (64) 



2.5.2. Matter Dominated Linearisation 

This is based on the radiation-baryon ratio, rj oc We keep every 0(ri) but in the 

almost-Friedmann-Lemaitre case of matter domination we then drop everything that is at 
least 0{r]v, rje, ev, rf,^^, v'^). We then have that matter dominated almost-Friedmann-Lemaitre 
means keeping terms that are at most : 



Matter dominated almost-Friedmann-Lemaitre ~ 0{e,v,ri). 



(65) 
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2.5.3. Expansion in Thompson Scattering Time 

Wc will introduce a pcrturbativc scheme in the Thompson scattering time, tc = {<TTne)~^, 
and will consider terms up to 0{t^) during the tight-coupling calculation such an expansion 
will be used to generate equations near to tight-coupling, the limiting case being when tc = 0. 
Additionally an equivalent scheme can be constructed in terms of the differential optical depth 
k' . This scheme is useful in the slow-decoupling era, i.e. in expansions where k' and (k')^ are 
sufRciently small to be ignored when compared to terms of order k. This approximation allows 
one to additionally consider the case when /t'e~" <C e~'^. 

2.5.4- Small and Large Scales 

We will find it convenient to introduce the notion of small and large scales. We will do this 
in two heuristically equivalent ways. The first scheme is based on the parameter ch, where the 
Hubble expansion is of order eh. and is used when considering situations outside the liubble 
flow; thus in the almost-Friedmann-Lemaitre small-scale case one would ignore all terms at 
least 0{e'jj,e'^,v'^,veH,V€,eeH). This scheme is useful since it can be used without a mode 
expansion, ft is ideal for making qualitative statements without the details which arise when 
introducing mode functions: specifically avoiding the complexity of mode-mode coupling in the 
small scale non-linear situation. By small-scale almost-Friedmann-Lemaitre we have in mind 
keeping only terms that are at most: 

Small-scale almost-Friedmann-Lemaitre « 0(eH, e, t)) . (66) 

A second and more precise scheme is based on the Hubble scale Xh = as/kH defined near the 
time of emission (E), allowing one to use k/kn > 1 and k/kn < 1 as characterizing large and 
small scales respectively. 

3. COVARIANT INTEGRAL SOLUTIONS 

This section has three aims: (i) Reproducing the integral solution of the free-streaming mode 
equations and modifying them in order to take into account Thomson scattering, using the CGI 
variables [40, 69]. We carry out a time-like integration, instead of a null-cone integration cor- 
responding to the original Sachs- Wolfe paper [68] , restricting ourselves to scalar perturbations 
with adiabatic modes only and assuming for the most part a, K = almost-Friedmann-Lemaitre 
background universe, (ii) Showing in the CGI formulation how fluctuations at last scattering 
time result in measurable cosmic background radiation anisotropies. (iii) Demonstrating how 
the solution can be related to standard formalisms by choosing specific frames; in particular 
we consider the Newtonian frame based on a shear free congruence. 

The basic equation we are concerned with in this section is the integrated Boltzmann equation 
(4). In covariant form it is given by 

f{x, e) + e"D„T(a;, e) + B{x, e) ~ C[x, e] , (67) 

where the gravitational source term, B, and Thompson source term, C, for damping by Thomp- 
son scattering are respectively given by: 

fi= -iD''T„ + (D„lnT + ^„)e'» + a«be°e^ C[x,e] ^ k{e''v^ - t) , (68) 
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and in the almost-Priedmann-Lemaitre situation in mind^^, the Einstein field equations give: 

3p-i(divi?)„«D,lnp„, aab--2[{p;,^Eab)-cn^\{pZ^Hai>)]. (69) 

The mode expanded form of this equation for the flat case {K = 0) can be written in the 
compact form as follows: 



r'e + k 



Te+i - T^-i 



k't^ ~ Se 



(70) 



where 



Sb = - [aBoSm + (aBi + k'vb)^« + aBafe] • (71) 

This combines (11-13) in a single equation for Te{ri,k) (see [18, 59]), where r = X]^r^Gf[Q] is 
written in terms of conformal time r] {dt = adr]), rather than proper time t. It is valid for all 
£ > 0, with To = a solution as required, consistent with the definitions we introduced above 
{Bo cancels the dipole term on the right in this case). 

In what follows, we will deal with the integral solutions to (70) given the source terms (71). 
The paradigm is to match an (almost-Priedmann-Lemaitre) era of free-streaming to one of 
tight-coupling. We will construct the homogeneous solution (without gravity or scattering) 
first, then include the gravitational effects to construct the free-streaming solution ( i.e. after 
decoupling to the present day) and finally include Thompson scattering to find the integral 
solution including scattering (which can be used during slow decoupling or to include effects 
of reionization) . Diffusion damping is included in this full solution with Thompson scattering, 
which in general has to be solved numerically, however it is helpful to introduce various ana- 
lytical approximations for the different stages described by the solution; this will be done later, 
where the visibility function approximation is used and the damping scale derived. Additional 
effects, such as the anisotropic correction and polarization correction, have to be dealt with 
separately. 

3.1. Integral Solutions (Flat Almost-Friedmann-Lemaitre Case) 
Here we wish to find the general solution to (70) without collision terms, i.e. with k' = 0, 
integrating along time- like curves using conformal variables. In order to do this we first find 
the solution to the homogeneous version of the above equation (i.e. for no gravity and no 
scattering), second, an integral solution of the inhomogeneous equation with gravity taken into 
account from the homogeneous solution, and third, the general solution for free-streaming. In 
the next sections wc consider the effect of Thompson scattering (k ^ 0), and the transition 
from tight coupling to free-streaming. 

The approach here is similar to the Seljak-Zaldariagga treatment [69, 83], however, they 
have taken the Sachs- Wolfe like formulation of the integrated Boltzmann equations, which is 
an integration down the null cone, and integrated out the angular dependence over Lcgendre 
polynomials (angular averaging) in order to construct the mode coefficients; then the conformal 
radial distance, X) is written in terms of the conformal time r], leading to an integral solution 
dependent only on the conformal time. Thus, formally they have carried out a null-integration. 
By contrast, what is carried out here is in effect an integration of the integrated Boltzmann 
equations down the matter world-lines, thus this is a time- like integration, onto an initial surface 
('last scattering'). The corresponding initial data near our past world line on that surface can 

^^One can compare this to the formulation of Durrer [16, 17] (eqns 3.5 and 18). To sec how this is linked 
to the Bardccn potentials ^ and we can use = ^D^„D5)(* — [17]. Notice that Durrer's integral 
solutions take on the Sachs- Wolfe form and can be compared with the treatments in [13, 61], while ours follow 
the form of [81, 69]. 
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then be related to data on the intersection of the past null-cone with the surface by means of 
a suitable homogeneity assumption. The time-like nature of the integration is often not made 
particularly clear in the literature, but solutions of the generic multipole divergence equations 
(from which the almost-Priedmann-Lemaitre mode hierarchy of temperature anisotropies are 
derived) are usually based on time-like integrations in the relativistic kinetic theory [59] . 

There is no effective difference between the time-like and null-integrations at linear order. 
This is because, at linear order, one needs to only integrate up one null geodesic, and then 
decompose the resulting temperature into its various moments. This is equivalent to making 
the multipole decomposition first and integrating these up one time-like world-line. This equiv- 
alence does not hold in the exact case, so when trying to include the effects of non-linearity, 
assuming such an equivalence could lead to misleading results. 

3.1.1. Finding the Homogeneous Solution 

The ^ = 0, 1 and 2 multipole divergence equations and hence mode form of the integrated 
Boltzmann equations are exceptional, given that tq = 0. The point of the integral solutions is 
to cast the exceptional equations ^ = 0, 1, 2 into a form that allows analj^tic investigation. We 
now find the covariant homogeneous solutions. 

Consider the homogeneous equation (valid for ^ > 1) : 



, (72) 



(^ + lf (0) ^(0) 



_ (2^+1) (2^ + 3)^+1 
for the background K = Q casc^^ without damping. The functions 

Tf\k,r,) = {2t+l)l3j^u{k7i) (73) 
are solutions of (86), if the coefficients obey the recursion relations [18]: 

(£ + = {21 + l)/3^+i, /3^(2^ - 1) = £/3,_i . (74) 

This can be shown by multiplying (72) through by = l3f (2£+ and comparing the result 
with the recursion relation for spherical Bessel functions. If the function jo{kri) satisfies the 
equation (72) for ^ = (there are of course no terms with i <1) then the rest of the equations 
(^ > 1) will be satisfied because of the recursion relations : 

-{2i + l){am)' ~ fc[(£ + l){at+iTt+i) - l{at-in-i)] . (75) 

The freedom in Piik) occurs in Poik) and (3i{k). Given that the Bessel function is finite at 
the origin: je{0) = 5io, these can be chosen to satisfy fHo = Pi = +1, the rest are generated 
through the recursion relations on Pi and then determine the solution r^""* (fc, rj). The arbitrarily 
specifiable initial data is later fixed by introducing an integral solution ((79) below) containing 
arbitrary functions Ca(?7) (see (80)) which are determined by the Einstein field equations 
through Bi{ri). 

The corresponding mode functions are 

oc 

r(o)(x,e) =Y,pj\2£+l)u{kii)0^'QA^,,r , (76) 
e=i 

(cf. (7)) and the corresponding multipole coefiicients can then be found: 

TA, U c^Pe'{2i+l) j> {kv)QA, U ■ (77) 



*Or for the open or closed cases by following [39] [35]. 
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Notice that this differs by a factor of from Wilson [81] since we are using plain mode 
functions instead of plane waves, although these can be easily related. Note there is no explicit 
mode mixing in this approximation, but such mixing is implicitly determined by the recursion 
relations (74). This shows that we should be careful with any truncation procedure we propose 
(see Appendix E). This procedure can be easily extended to the open case using the recursion 
relations for the open mode functions (see Appendix K). 

3.1.2. Construction oj the Integral Solution 

Given that we have the solution Tf^\k,ri) to the homogeneous equation of the form (86), 
now consider the equation with given gravitational source terms, but still without damping: 

(^+1)2 



r'e + k 



{2e+l){2£ + 3) 



n+i — n-i 



= - [aBiSn + aB2 6^2] ■ 



(78) 



What is important to notice here, is that this equation is valid for £ > 1, not ^ > as in (70); 
indeed tq = 0. We need to find a particular solution to this equation. 

We proceed as follows. Consider the ansatz in terms of using 6r] = rj — rj' , along with the 
Liebnitz rule for differentiation of integrals: 



rfiv) 



dr]'Ae{r],r,') -^{rj) = drj' —Ae{ri,r]') + Ai{r],ri) . (79) 



Now we define the kernel, A^, as in [81, 83]: 

.(0) 



(0) 



dr," 



{Sr,) 



It can then be shown from (79) and (80) that : 

(^ + 1)2 



(2£+l)(2£ + 3) 



Sc, 



(80) 



(81) 



where 



Sc = Co(r?)rf (0) + Ci(j?)rf (0) + C2{v)t^'^" {0) , 



(82) 



given Tf{ri) as in (79-80). Thus (78) is satisfied by our ansatz provided the coefficients 
Co{r]), Ci{r]), and C2(?7) in the integral solution are found in terms of the CGI variables 60(^7), 
Si (77), and ^2(77) determined by the Einstein field equations. This will be considered next, 
when we put the parts of the solution together to obtain (90). 

3.1.3. Inclusion of Damping 

Here we extend the previous solution (78), where the relationship between the coefficients 
in the integrated Boltzmann equations (D.l) can be read off from (82), including damping 
through k' . We notice that if rf'^ is a solution to (86), then 



(83) 



will be a solution to 



{£ + lf 



(2£+l)(2^ + 3) 



I / * 



. 



(84) 
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Similarly we find that for the integral solution (79), of (78), the expressions 

rf{ri)= r dri'e-<^'^Ai{ri,r,') or r^^ (r^) = e-^^^^ [' dn'Ae{v,v') , (85) 
Jo Jo 

will be particular solutions to (70), given the correct choice of Co, Ci and C2. Hence we 
can modify the solutions of the previous section to include Thompson scattering by simply 
including the damping terms as in these equations. Hence the extended equations include the 
special case of free streaming, when for some interval of time k = 0; thus they can extend all 
the way from late tight coupling to the present day, if we include a suitably time-dependent 
scattering coefficient k. 

In more detail: we have that T^^\r]t,) = 0. This is simply due to the fact that during tight- 
coupling there are no higher moments, just the monopole. Here we assume that free-streaming 
begins after some ?7* . The slow decoupling solution will modify this assumption. As before we 
have (now using Sr]* = ?? — ??*) 

(r;) = {2£ + l)(3^'je{k6v*) ^ t;^(7?) ^ / dv'Ae{r^,rj') , (86) 

Jo 



where the initial data is now given by Ci{r]' — r]^,)\ notice that we do not introduce an additional 
77* as we will be using the solution of t^^^ already including the initial conditions Once 
again we have the integral solution t[ (r/) integrated from to 6ri* such that the anisotropics 
are now determined by t[ {r}o)- Evaluating the integral from to A?7* (for A?7 = r]o — r]'), we 
find 

rvo r- 

Jv* 

+ Ci(77'-77*)rf)'(A,?)+C2(V-»?*)Tr"(A»7)] , (87) 
where the damping term now enters explicitly. 

3.1.4- The Complete Solution 

We can now construct the general solution to (70) with 7^ by putting the previous 
results together. The homogeneous seed solution T^^\r]) is given by (73). The particular 
integral solution rf (jj) is given, in terms of Sr] = r] — r]' by (85). The general solution is then 
given by 

r,(,7)=e-«Wrf)(r?)+r;^(77). (88) 

Substituting this into the general equation (70) and using the radial eigenfunctions evaluated 
at zero (in particular je{0) = Seo along with the recursion relation) gives 

CoSm + Ci p^^kSei + C2kHl3^^ (^^5^ + ^<5^o) = - [(aBi + K'vB)5n + ^^2^] , (89) 

relating the functions determining the solution to the time-dependent coefficients in the equa- 
tion. From (89) the functions in the integral solution are found in terms of the dynamical CGI 
variables: 

-Co{ri) - +\aB2. -C^{^)c^+\{aB^ + n'vB) , - C^iv) - +^aB2 . (90) 



^®We could have used rj'^-'(<5??) Ti'^\Sr] + r}^) along with the original homogeneous solution unchanged. 
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For scalar perturbations, the term C2 is ofFcctivcly the coefBcient of the shear; the term Ci is 
the coefficient of the gradient of the temperature and the acceleration and the coefficient Co 
also contributes to the shear (see (G.19)). Note that these are CGI with respect to 

These functions are all evaluated at a time r], which takes all values from rid, the time of 
decoupling, to 770, the time of observation. In the nuU-conc formulation of the problem this 
input of new information corresponds to the way the null geodesies from the point of emission 
to the observer keep crossing new matter and hence encounter new information. Because we 
are integrating on a time-like curve, this information is represented here as varying with time 
along that curve; and in some simple circumstances, the values at later times are determined 
fully by the values at earlier times (as happens, for example, in the original Sachs- Wolfe case: 
K = 0, p = 0, and only growing scalar modes are considered). 

3.2. Integration by Parts 

In order to deal easily with the initial data it is now useful to write the general solution for 

K ~ in terms of the present time, 770, and the initial time, 77*, by integrating with respect 
to conformal time and defining Ary* — ■ We could fix the conformal time by setting 
r]o = here. In order to recover the results of [10, 83] one would take 77* 0, however, we 
would like to recover the results as close to [39] and so retain their conventions where possible. 

Notice that from dr^^^ /dt = 0, we have r'^^ (x(?7), e(?7)) ~ T(°)(a;(77*), e(?7,)), and ^'^(0) = Seo 
(we have chosen the solution to be finite at origin). We choose the initial conditions Te,{rii,) = 
[81] and ro(»7) = and using the parameter freedom in the homogeneous solution, set Po{k) = 
+ 1 and Pi{k) = +1. The homogeneous solution is now fixed as in (86), for Sr]* = rj — rj^, and 
this in turn sets the integral solution to (85) [81]: 

rl'^iv)^ e-'^^^^Mn,v')dv' , (91) 
Jo 

where we still have the freedom of setting the initial data for the integral solution from the 
C/(r7)'s which are fixed by the Einstein field equations. Putting this all together we find 

reivo) = r^"\v*) + rfirio) = r;^(77o) = / e-'^ Ai{r}o,n')dr}' . (92) 

Jo 

On changing the integration to from rj., to rjQ in (87), integrating by parts, and using the 
initial conditions (once again T^^°'(r7*) = 0) we find : 

T,(r7o) ^ [C'M - Ci(7?.)] e-'^rf ^(770) - C2(7?.)e-'*rf ^'(770), 

+ / d77'e-«[Co(77')-Ci(77') + C^(r/')]Tr(% + r?*-77'), 

+ / dv'in'e-^) [Ci (77') - 2C2(7?')] ^ (770 + 77* - 77'), • • • 

Jv* 

rVo 

+ / dr,'{{K'f-K")e-^C2{v')T^°\m+V*-v')- (93) 
Jv 

The initial data for the solution Tf^\k, rf) is the set of constants Ci{ri^) which are determined 
by BA{k) = {Bo{k) , Bi{k) , B2{k)}; these must be matched to the initial distribution function 
on an appropriate initial surface S (for example, the 'surface of last scattering' which can be 



^''The relationship between the conformal time 17 and the radial distance x is dx = ~'^'7 so x = V ~ V which 
follows for the homogeneity and isotropy in the background. 
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covariantly and gauge invariantly defined). This then determines the solution up to the present 
day (and after). We are free to chose any Q's as long as they solve the Helmholtz equation 
in the background. The choice of Q then explicitly determines for example we are free 

to choose Q to be the spherical or plane- wave basis. In practice we naturally use two sets of 
mode fimctions Ge[Q], matching those for the nuU-conc (given in a spherical basis) to those 
in some initial surface (given in terms of a plane- wave basis) . The matching of these two sets 
of harmonics is then given by the relations usually written into the construction of the mode 
coefficients (see (76)). This matching is based on mode functions Gi[Q] in the Robertson- 
Walker background, which is acceptable because of the homogeneity assumption. By using 
Ge[Q] we do not actually need the explicit form of the Q's. 

Equation (93) shows (r.h.s. of the first line) how major parts of the cosmic background 
radiation anisotropy arc determined directly from the set of initial conditions (at last scattering, 
for the freely propagating radiation). The integrated effect arises through the coefficients Ci{r]) 
as integrated down time-like geodesies in the remaining terms on the right hand side. In general 
there is a non-linear coupling through the field equations between the matter, the radiation and 
the acceleration and shear terms that arise in the integrated part. The situation is much simpler 
when this back-reaction can be neglected; for this reason it is convenient to consider the case of 
matter domination, during which the radiation can be considered as a test-field propagating on 
the background determined by the matter content. However we can also consider the general 
set of linearised field equations (see Appendix G.l) and the coupling to the radiation via the 
source terms, first the gravitational source, B, and second, the scattering source, C[r], (68), in 
(67). In the following sections we look at the various approximations that can be applied at 
different epochs. 



4. FREE-STREAMING 

Using the integral solution (93) we construct the almost-Friedmann-Lemaitre free-streaming 
projection of the initial conditions near last scattering to here and now (the determination 
of these initial conditions is demonstrated in latter sections) and the integrated secondary 

contributions arising during the period after last-scattering until now (we have dropped the 
baryon relative velocity effect using the instantaneous decoupling assumption): 



{2e + 1) ' 



■aBo 



1 



{aB^y 



5 1 
3fc2' 



(062)" j^(fcAr,) 



(94) 



Where we used as final conditions : 



[aSi]'(77o) = [aS2](7?o)=0. (95) 

The first term on the right, the Bi term, will generate the acoustic primary effect on the 
anisotropies, the second term is the Doppler contribution due to the radiation dipole (the 

baryon velocity contribution which would arise through Ci (90)), the third and fourth terms 
give the effect of any shear, near last scattering (through the initial conditions of S2). The 
remaining terms represent the integrated Sachs- Wolfe effect. 

The above equation will be modified in the following section to include slow decoupling, but 
first we demonstrate how to recover the basic Sachs- Wolfe effect. 
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4.1. The Almost-Friedmann-Lemaitre Sachs- Wolfe Effect 

We now find the solntions corresponding to the matter dominated, free-streaming era, with 
adiabatic modes only, using the Newtonian frame treatment (see section 7.1.2). Using the field 
equations from [59], or from the Appendix G.l and 59, we can find the source terms aBi{k,t) 
for the free-streaming projection : eqns. G.24 through G.26 in Appendix G.2. This applies to 
the case of instant decoupling. 

During matter domination the dipole is negligible, so we ignore it. The shear contribution 
is small on large scales, hence we can also ignore it. 

On substituting these equations into the flat almost-Friedmann-Lemaitrc integral solution 
(93) with K = in the source terms (eqns. G.24 through G.26) we can find the free-streaming 
almost-Priedmann-Lemaitre solutions for the temperature anisotropics (94). Using ($/9~^)' 

we find the CGI kinetic theory equivalent of the Sachs- Wolfe formula for cosmic background 
radiation anisotropics in terms of matter inhomogcneitics at last scattering at various wave- 
lengths, together with an integral term. In various cases (see [41] for references to more general 
treatments), in particular the matter- dominated spatially flat solutions with only growing scalar 
modes, the integral terms vanish and we obtain: 

This gives the (approximate) projection of the large scale potential inhomogeneities at last 
scattering onto the sky today; the mean-squares jr^p can then constructed, using the results 
from Paper I [18]. The effect arise from the terms DalnT « ^Dolnp « p~^D^Eab, having 
used the adiabatic assumption. 

This recovers the standard Sachs- Wolfe result [68, 13] on large scales, where the potential 
fluctuations are just due to primordial initial inhomogeneities that are unchanged by interven- 
ing physics. 

We discuss the origin of these fluctuations in later sections - they are given by solving 
these equations before decoupling, which for example implies the existence of acoustic oscilla- 
tions. These potential fluctuations are what seed structure formation through the production 
of matter perturbations undergoing gravitational collapse beneath the Jeans scale. The mat- 
ter perturbations effectively decouple near matter-radiation equality, making the large scale 
temperature anisotropics the key link between the radiation anisotropics now to the potential 
fluctuations then (near radiation decoupling), and so to the matter power spectrum both on 
large and small scales today. 

Notice that these equations will hold for any choice of 4-velocity that is close to the matter 
4-velocity, i.e. there is still frame-freedom associated with this freedom of choice. As has been 
remarked various times, there are several possible physical choices for this 4-velocity (which will 
all agree at late times); the interpretation of the physical meaning of the cosmic background 
radiation anisotropy sources will change depending on this choice. The important difference 
about the derivation of the Sachs- Wolfe effect here as opposed to other treatments is that (i) 
this result is found in the total matter frame, (ii) the integration is explicitly time-like, so it is 
not treated mathematically as a projection along null rays but rather as the evolution of the 
anisotropics of radiation in a small comoving box, as explained in the introduction. Thus the 
initial data here is not at the intersection of the past light cone and the last scattering surface, 
but rather at the intersection of the world line of the observer and the last scattering surface. 

This analysis can be compared to primary anisotropy source term of the gauge-invariant 
treatment used in [40]. The subtle difference between the Bardeen variable gauge-invariant 
approach and the CGI approach used here is that the Doppler source, which in their case 

^''In fact one has that for an Eistein-de Sitter background : [^p~^]{ri*) = —{3H^Qo)~^ D*[k'^^A{k,0)] and 
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arises through Bq, now enters through the integral terms only; there is no direct Doppler 
contribution at last scattering from the first term in the integrated solution. 



5. SLOW DECOUPLING 



To deal with slow-decoupling, we return to the general damped solution, (93), and introduce 
the slow decoupling approximation 



«"<1, («') <1. 
to find, on substituting in from the coefficient relations (90): 



(97) 



{2£ + 2) 



1 n° 1 

+-(K'e-'')ws(?7*)i^(fcA77*) + j drj'iv'gK' + K"vB)e-'')-je{kAri) 
rm r 1 1 



je{kAr]) . 



(98) 



We see that damping effects are controlled by e~'^ and re'e"". A further approximation would 
be to take e~'^ K'e~'^, so that we need only consider the free-streaming like solutions, which 
we then convolve with the damping factor, defined as a combination of the visibility function 
and the diffusion damping envelope - this is done later using the damping envelope as derived 
from the dispersion relations in section (5.2). Later we will explicitly recover these equations 
in the Newtonian frame of section (2.4.1). 

5.1. Silk Damping 

Diffusion damping will occur and introduce a damping scale, the Silk scale, giving a cut-off 
in the matter perturbations, and there will be a corresponding diffusion damping effect in the 
photons. This is naturally included in our general damped solutions in terms of the exponential 
envelope implied by the equations, which can be demonstrated heuristically. 

The cut-off arising through photon diffusion occurs when the term involving k' in (70) 
dominates the other terms; that is, when for any £, k is large enough that 



k 



(99) 



the approximation assuming that the damped modes are roughly of the same magnitude 
(independent of £ when this condition is satisfied). This then implies an exponential decay in 
the relevant modes: 



4 

k -C -k' 

o 



W -K Te 



Teiv) ~ exp(-K'77)r^(0) . 



(100) 



Thus small scales will be heavily damped by this process and long wavelengths unaffected, 
leading to a wavelength-dependent damping envelope. The resulting cut-off in perturbation 
amplitude at a critical wavelength at last scattering will result in a corresponding cut-off in 
cosmic background radiation anisotropy amplitudes observed at a critical angular scale. 

A more detailed examination undertaken later will show the explicit wavelength dependence 
of this cut-off effect. 
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5.2. The Visibility Function 

An alternative approach to the slow decoupling solution (98), is to argue that the dominant 
contribution during slow-decoupling arises from the visibility function defined by V{k,r]) « 
K'e~'^ as convolved with the free-streaming integral solution. The visibility function gives the 
probability of a photon last scattering during a small time interval drj. From Hu & Sugiyama 
[39] it is useful to define the damping factor (now including diffusion damping, which will be 
derived from the coupled baryon-photon equations in (169) ): 

rVo 

V{r]o,k)= rf?7V(r/,A;)e-('=/'=^)' «e-('=/'=°'\ (101) 

The visibility function will model the changing ionization fraction, this does not include the 
diffusion damping, which is added in by hand above, through the damping scale. It should be 
realized here that the Gaussian diffusion damping is naturally included in the original almost- 
Priedmann-Lemaitre integral solution (93). However, given that we will use solutions that 
are first-order (in the scattering time) and then recover an explicit dispersion relation for the 
damping scale, at second-order in the scattering time, it is convenient to modify the damping 
factors such that they are re- written in terms of the visibility function 101. The second order 
damping scale of the form used here is explicitly derived in section 6.2.5. 

Now, we can modify the free-streaming projection by including the damping factor V and 
the baryon velocity effect (which must be put in from (93)). In this approximation, we can 
effectively drop the last two lines of (93) except for the initial baryon velocity contribution, to 
obtain: 

re{vo)Pi ^ [c'^^rj.) - cM]v{rjo,k)MkArj,) 



(2^+1) 

+ fcC2(r?* )!>(%,/=) 



-je-i{kAr]^) - ^\ ji+i{kAri*) 



(2£+l)-" ' ' {2£ + iy 

+ / d77Ve-('=/'=-) {Co{v)-C[{r,)+C'^{r,)}MkArj), (102) 

where the coefficients Ci , C-z and C3 are given by (90). This has the effect of taking the previous 
more general solution (93) and specializing it to the most important regime as far as decoupling 
is concerned, thus giving a major improvement on the sharp decoupling approximation, while 
avoiding the complications of the complete integral solution given above. A similar correction 
is made using the visibility function in the integrated part of the solution, in order to best 
deal with a changing ionization fraction, given that we will once again only be using almost- 
Priedmann-Lemaitre solutions, that are either first order or zero-th order in the scattering 
times (discussed below). 

5.3. Slow Decoupling in the Conformal Newtonian Frame 

Here we cast the above derived solutions (94,98, 102) based on the integral solution (93) into 

the CGI Newtonian frame (based on the shcar-frcc frame described in section 2.4.1) for the 
case of scalar perturbations, in terms of the Bardeen like scalar potentials and 

The vanishing shear condition aab ~ implies that C'o('?) ~ and C2{r]) ~ 0, hence we 
can use these conditions directly to find the slow- decoupling anisotropy solution for an almost- 
Friedmann-Lemaitre model in the Newtonian frame: 

0, -C,{7j^)V{rjo,k)M{kA7j^)- r d7jVe-^>'/'^-^'c[{v)M{kAr^). (103) 



I?. 
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Using the results from section (2.4.1) and the almost-Fricdmann-Lcmaitre relations (given in 
Appendix G), it can be shown that the key quantity of interest i?a ~ Do In + ^a^A, in the 
case of scalar perturbations, obeys the following relation in the Newtonian frame. 

£>( Jb) « D(„Db>($^ - - 2ifD(„Db)$A - ^D<„D6)(D^fe) (104) 



Now we need to find the mode coefficients for B in terms of the quantities defined in the 
Newtonian frame. Writing 

D„lnT « ^(5TQ„, D„$a « f^^Qa, and « f^^Q^ , (105) 

we find that 



(106) 



On mode expanding (104) '■^ and transforming to the conformal time derivative we obtain: 



{aBiY « -a^HBi + fc($'^ - ^'h) - 2Hak^A + \k'^fi- 
Now using (90) we find that: 

-C'M-in'vB)' ^+\{aBr)'. 
We can now put this all together, first, from (90) and (106) to find: 

-C^{ri)^K'vB + {5f + ^A). 
and second from (108), (109) and (107) to find: 

-C[{ri) « {k'vb)' + ({^'a - ^h) - aH{5f + ^a) - 2aH^A + |fcfi) 
Substituting these results into the integral solution (103) we obtain: 



(107) 



(108) 



(109) 



(110) 



{2£+l) 



{6T + $a) + k'vb {v*,k)V{r]o, A;)j>(A;Ar?*) 



+ 



drjVe 



-{k/kof 



{k'v'b + k"vb) + lkfA jeikAri) 



+ J dryVe-('=/'=") {{^'a - ^h) - aH{5f + 3<^a)\ jeikAr^) . (Ill) 

Here the second order terms (both in terms of the scattering time and in the almost-Priedmann- 

Lemaitre sense) have been dropped. Wc can then pull out the canonical solution when wc ignore 
the Dopplcr contribution, the initial baryon relative velocity at last scattering (it is tightly 
coupled to the radiation velocity and is thus small already). We also ignore the intermediate 
scale integrated effect which contributes to the early-integrated Sachs Wolfe eflfect. The result 
is: 



{21 + 1) 



5T + ^, 



T^im, k)ji{kAr]*) 



Vo 



dT^Ve-^^/^^^^ [^'a - $y j£(A:A?7). (112) 



8We will use D<„D(,)iyi « -^A^{a^b)Q, D"Tc « +|riQ and D(„B(,) 



t(Bi+2m)D(„Di,>Q. 
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This completes the recovery of the standard integral solution results using the 1+3 CGI 
approach at linear order - we have notationally suppressed the fc-dependence of the temperature 
anisotropy (t^(??o) = Ti{rio, k)). It corroborates the standard anisotropy derivations based on 
a 3+1 hypersurface fohation, which uses the Bardeen formalism in the conformal Newtonian 
gauge. 

6. LATE TIGHT-COUPLING 

Here we extend the Thomson scattering analysis of the previous sections to include a simple 
model of late-tight coupling and hence of fast decoupling. 

We aim to reproduce in covariant form the Peebles and Yu near-tight coupling [66] and Hu 

and Sugiyama tight- coupling approximation [40] treatments, valid for the period of late tight- 
coupling, up to and including decoupling Remember that we are ignoring the anisotropic 
and polarization effects as these can be corrected at the level of the damping scale. 

6.1. Integrated Boltzmann Equation: Near-tight Coupling 
Here the almost-Priedmann-Lemaitre integrated Boltzmann equations is used to construct 

a set of multipole divergence equations that describe the radiation near tight-coupling. These 
are the intermediate scale equations, valid in the tight-coupling era. This is done by carrying 
out a CGI version of perturbation theory in terms of the scattering time. 

6.1.1. The Scattering- strength Expansion 

The solutions we have considered so far are linearised through a small-parameter expansion 
in terms of the anisotropy parameter r. The basic idea now, following the method of Peebles 
and Yu [66], is that additionally a second expansion is constructed in terms of the collision 
parameter tc = {(7x^1,;)^^, without truncating the Boltzmann hierarchy at the order of the 
calculation, and thus avoiding the problems inherent in exact truncation [29] (see Appendix 
E). We thus find the evolution equations for the energy density, momentum flux, and the 
anisotropic flux of the radiation close to tight coupling. 

Consider the almost-Friedmann-Lemaitre integrated Boltzmann equations (4) and (67) for 
isotropic (in the baryon frame) Thompson scattering (68) [59] ; this is inverted to find : 

r(a;', e°) = w^e„ -tc[B + f + e"D„r] . (113) 

We now systematically approximate (113) in terms of the smallness parameter tc- The right- 
hand side (the scattering term) is used to find the zero-th order collision-time correction to the 
total bolometric temperature, with corresponding temperature anisotropy given by 

r(o)(a;\e'')«z;ge„ . (114) 

The equation is now perturbed about the zero order velocity perturbation and one can then 
recover the first and second order corrections in tc to the zero-th order temperature anisotropics, 
to find, r(i) and r(2) respectively, where the n-th. order correction is denoted by r(„). We obtain 
an almost-Priedmann-Lemaitre perturbative expansion in tc ■ 

T(n) « V%ea - tc [B + f(„_i) + e"DaT(„_i)] . (115) 

^^Here we are explicitly making a distinction between the treatment [40] (what we call tight- coupling approx- 
imation) and that in [66] (what we call near-tight coupling. By tight-coupling approximation we mean that 
is sufficiently small that it can be ignored (inducing a contribution of the order of magnitude (say) of at 
least 10~^) when multiplying quajitities of linear order such as the shear) The near-tight coupling includes the 
radiation quadrupole in the case of isotropic Thompson scattering (in the matter frame). 
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The tight-coupling Hmit is recovered when tc = 0. This treatment is then a consistent (in the 
sense of the truncation conditions described in Appendix E) near-to-tight-couphng treatment 
in almost-Priedmann-Lemaitre universes. (The first, in n, three temperature anisotropies, are 
given in Appendix C). 

6.1.2. Solid Angle Integration 

Now the temperature anisotropy is integrated over the soHd sphere to ensure the condition 
that there is no contribution to the bolometric average T((,)(a;*), 

I T{x\e"-)dn = . (116) 

It should be clear why the second order correction to the temperature anisotropy is needed 
even though we intend to keep the expansion only to first order in t^, the integrations over 
term the e"'Va will vanish. Now by integrating T(2) (C.6) over the solid angle and using (116) 
and orthogonality of O^* the gradient of the radiation flux is found: 

D„r» ~ T)av% - ac [(DaV%) - (D^r") + \B\\npR) + {T>aV%) + (D„A«)] . (117) 

By taking spatial gradients of the radiation flux (121) we flnd on comparing with (117), that 
in order for there to be no contributions to scalars, 

(divi;)^ « (divr) . (118) 

The above equation (117) then becomes 

(divr) ~ (divt;)B - tc [iD^lnpfl + (D„i)^) + (divA)] . (119) 



6.1.3. The Transport Equations 
Finally the individual PSTF multipoles are recovered at a given order 



a; 
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0^,T(„)(a;%e«)dO. 



(120) 



Here is defined as before in Paper I [18]. The second order temperature multipoles arc now 
found from (120) and integrating T(2)(x, e) (C.5) after the combination of direction vectors has 
been replaced by PSTF tensors (one can use C.7-C.9) : 



(D^ InT + A^) + i)^ - iD''DcT' 



-ab 



^abc 



(o.bc) 



(D<"z;^>y + D<"(D^> InT + yl^>) + D<"4^ + tr"' 



>3 



(121) 

(122) 

(123) 
(124) 



where we have dropped terms of 0{il). These are the key results of this section. They are 
the appropriate transport equations for the late-tight coupling era, i.e. up to last scattering, 
and are essentially equivalent to the causal transport equations given by causal relativistic 
thermodynamics [60]. 
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What wo have shown hero is that if wo are interested in the behaviour of the photon-baryon 
systems to first order in the scattering time, a dissipative fluid approximation is sufficient 
to describe the radiation (cf. the papers by Israel and Stewart [45]), and will not result 
in an explicit truncation of the Boltzmann multipole hierarchy, rather it gives a systematic 
approximation scheme where wc can, to the appropriate accuracy, ignore the third order and 
higher terms. This is significant; one cannot merely drop the higher order moments and 
truncate to a fluid description, as the kinetic theory treatment fixes the transport equations. 
Here we have consistently decoupled the I < 3 multipole equations from the rest of the hierarchy 
and the consistency of this decoupling is maintained through (117) and (121 - 124). 

The solutions to these equations, which lead to acoustic oscillations during this period, will 
then affect the cosmic background radiation anisotropics by setting initial conditions for the 
free-streaming solution discussed in the previous section. We give a derivation of these results 
in the following section. 

6.2. Late Tight-coupling and the Oscillator Equation 

Here we derive the CGI equivalent of the analytic tight-coupling approximation used by Hu 
& Sugiyama [39, 40] . This approach uses the tight-coupling limit in order to get rid of the radi- 
ation quadrupole during late tight-coupling ^° and covariantly reproduces Hu and Sugiyama's 
conclusions about the cosmic background radiation anisotropy due to inhomogeneities, acous- 
tic and Doppler sources (what they call 'primary sources'). This gives the 'Sachs- Wolfe effect' 
due to the Newtonian potential near last scattering, but not the 'integrated Sachs- Wolfe effect' 
due to changing potentials after tight coupling (resulting from more complex matter models 
and/or spatial curvature). 

6.2.1. Near Tight- coupling Equations 

We start with the near-tight coupling equations (121, 122) and (124) except rewritten to 
first order and in terms of the optical depth so as to be closer to the notation of the better 
known treatments [40] : 

To ~ - k-^ [D„(lnT) + Aa+ v^] , (125) 

Tab =i -K 

T^' ~ y£>2. (127) 



(126) 



Here we have assumed that the collisions are dominated by Thompson scattering and is 

therefore isotropic in the matter frame. 

The relative velocity of the matter with respect to the preferred reference frame is v'j^ ~ 
u% — u°'. Rewriting (126) in terms of the shear of the baryon frame, we have Tab — —f^~^^{aU^y 
so the quadrupole is given entirely by the shear of the matter. 

Notice that iTab = pRTab ~ 0, as the case of matter domination. This condition is not 
sufficient to ensure that Tab can be ignored in equations when it appears on its own, even 
though the quadrupole is small. The key point, which was discussed in section 2.5, is that 
there are four principal linearisations: the almost-Friedmann-Lcmaitrc one at least O(e^), the 
almost-Fricdmann-Lemaitrc radiation isotropy one 0{er]), 0{rf), (implying the previous by 
the almost-Ehlers-Gercn-Sachs theorem), the non-relativistic assumption 0{vri), O(w^), 0{ve), 



^^li there are any anisotropic contributions such as anisotropic scattering (in the matter frame) or large 
shear (from gravitational waves) at that time this sort of approximation should be considered with care - such 
phenomena would break tight-coupling. 
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and the linearisation scheme based on the differential optical depth. Hence care must be taken 
when approximations are made to the equations. 

6.2.2. Tight coupling: Momentum Equations 

The tight-coupling calculation is now carried out, assuming (125-127) hold. Consider once 
again the radiation energy and momentum conservation equations (£ = and £ = I multipole 
divergence equations): 

(lnT) + ie ~ -^divr, (128) 

-fa ~ A + Da(lnr) - k{v!^ - Ta) + ^DV^c , (129) 



and the baryon energy and momentum conservation equations 

(lnp) +e =^ -D»z;f, (130) 

-V^ ~ +HV^ +Aa+R-'k{v^ -Ta) . (131) 

Here k is the optical depth, and the radiation-baryon ratio in the real universe is given by 
(using the enthalpy h = p + p) by 

Ri:c^) = ^ = EK+PE ^ 3^ =,Rc,HR. (132) 

hR{x) PR+PR ^Pr{x^) 

This is related to the speed of sound in the background via = (l/3(iio + !))• The matter 
momentum equations, (131), give 

v^^ra-^[v^+Hv^+Aa] . (133) 

Substituting (125) into (133) and retaining all terms up to linear order (in the relaxation time) 
we obtain 

t;f ~ r„ - Rk-^[fa + Aa + Ht"] + Oik-^) . (134) 

This is then substituted into (129) in order to remove the velocity terms, and with a little 
algebra, we find 

-fa c.ua + TT^D«lnT+ - ^ (^D^D,„.^, (135) 

where the last term has been written in terms of the matter shear. We can now consider 
the situation where k~^ becomes sufficiently small (but non-zero) so that the last term can be 
ignored. This is possible as the matter shear is already at least first order. We then find (c/ 
[40] eqn. (B2 b)): 

• ^ ^ -DailnT) -Aa . (136) 



{1 + R) " {1 + R) 



^^This is found to 0[1] by substituting the matter energy conservation equation (dust part of G.IO) into 
matter momentum equation (dust part of G.ll) all to 0[1]. 

^^By speed of sound we mean adiabatic sound speed : = -■ — — = i-?- — '^r^ — r = ^t- — 7^ — vtt = 

^ I- s PM+PR 3 ipM+PR) 3 {pM/PRJ + i 

for matter domination Cs 0. 



3 (3pM/4.pR)+l' 
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This is the momentum flux equation for the radiation and is a key result. It can be rewritten 
as 

[(1 + R)Tai + D„(lnT) ~ -(1 + R)Aa , (137) 
or on taking its divergence as 

[a(l + i?)(DVa)]- + a(D' InT) ~ -(1 + i?)a(D"Aa) . (138) 



6.2.3. Spatial Gradients and the Oscillator Equation on Small Scales 

The basis of this derivation is the 'smafl-scale' assumption which effectively means that on 
small enough scales we can ignore the expansion (see section 2.5). This is just the statement 
that the scale of inhomogeneity is less than the Hubble scale (66), so we drop all terms of 
O(eeff) (see section 2.5.4). 

Our aim is to recover the standard oscillator equation (the somcc equation for the acoustic 
oscillations) using the 1+3 CGI formalism. Note however that we still have the freedom to set 
the relative velocity in the small boost equations (which we will do in the next section). 

Taking the spatial gradient^^ of the radiation energy conservation equation, (128), we find 

-iD„(D,r^) :^ (D„ InT) + ^D^e + H{B, InT + A,) , (139) 
and taking the divergence of the resulting equation above gives 

-^(D2(Dct'=) ~ {D^lnTY +2H{D^\nT) + \D'^& + H{divA) , (140) 
and this can be written as 

-ha'^B'^i'D^T^)) ~ (a^D^lnT)- + ^(a^D^G) + if (a^div A) . (141) 
This is analogous to equation (B3) in [40]. Then using (138) we find 

[(1 + i?)aDe(D„r")]' + H{1 + i?)aDe(D''Ta) + 00^(0^ InT) ~ -(1 + i?)aDe(div A). (142) 
Substituting (139) into (142) and using the fact that HB'^Aa « 0{eHe), we obtain 

-3[a(l + i?)(aDc InT)' ]' + a^Bc{D^ InT) 

~ [a(l + R){aBc&)] - a(l + i?)Dc(div A) . (143) 
Finally, transforming to conformal time, dt = adr] gives 

3[(1 + R){aBc InT)']' - a'^B^{B"-{aBa InT)) 

~ -[o(l + R)(aBce)]' + a2(l + i?)Dc(div A) . (144) 

On using the small-scale linearisation scheme described in section (66) we find, on dividing 
through by 3(1 + i?), the oscillator equation: 

(aDe InT)" « +-^-^D,(D2 InT) - a'{B,ey + yD,(DM„) . (145) 

23 Using the identity (Da/)' ~ D^/ - HDa/ + A^/ we find that (Da InT)' ~ DahiT - if(Da lnr + A^) from 
the almost- Priedmann-Lemaitre £ = multipole divergence equations and H = ^ [59] . 
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This is the covariant harmonic-oscillator equation which describes the acoustic modes We 
can compare it to the usual gauge invariant result found in the Newtonian gauge by transform- 
ing to the shear-free frame = where D^^n;,) = 0. 

We will investigate this equation in more detail in the next section and relate our results to 
those in the standard literature (which are expressed in the Newtonian gauge). 

6.2.4- The Newtonian Frame Oscillator Equation 

In this section we recover the harmonic oscillator equation of Hu and Sugiyama in the 
Newtonian gauge [39] from the CGI formalism. The difference between this and the previous 
section is that here we apply the small scale approximation at the very end of the calculation. 

Using the Newtonian frame choice ~ + v^, D^arif,) = 0) equations (136) and (15) 
become 

^ ;^a + T^^a(lnT)«-ia«-Da$A, (146) 



" 1 + R 1 + R 
{Da InT) + H{ba InT + i„) + ^D^Q « -h)a{D,f^) . (147) 

These can be re-written and put into the following form by using the almost-Friedmann- 
Lemaitre Einstein field equations, together with the transformation relations given in Appendix 
G: 

{Dan+H{J^af<^) « -^(D^f'') - ^{b^inT) - {D^^A) , (148) 
D^,D,)i\nT) « -D^„Dfc)$^ - ^B^,B,){D,f-) . (149) 

Taking the time derivative of (149) and substituting into equation (148) after first taking 
PSTF derivatives we obtain the full equation for D^^D^^ InT: 

(D(„D6) InT)" + £r(D<„D(,>$^ + 2D(„Db> InT) + H [(D<„D6)$^) + 2(D<„D6> InT)] 

« -{D^,B,^^h) -2HB^,B,^^h - 2if(D<,Db)$) + _2_D^„Dfc)D2lnT 



- B^aB^D^'S^A 



^ 3H 



(l + R) 



(D<„Db)(lnry + D(,Db)$a) . (150) 



On dropping all terms 0{eeH), we once again obtain the 1-1-3 covariant form of the small scale 
Newtonian frame oscillator equation (without using a mode expansion): 

(D<„D,)lnT)- + ^-A_D<„D6)(lnT) + ^D<„D,)D2lnT 

« -(D<„Db>$^^) + ^^^^B^.B.^^H - D(„D(,)D2$^ . (151) 

The techniques used to derive the above equation become useful latter when dealing with the 
non-linear terms as they avoid the complication of mode-mode couplings when understanding 
the qualitative features of various effects [59] . 



^"^This can also be obtained from a two-fluid CGI description [13], as well as from the imperfect-fluid descrip- 
tion [60] - the point here is that we have derived it from a self consistent kinetic theory approach, listing along 
the way the necessary physical approximation required to reduce it to the standard acoustic oscillator form. 
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Upon using the mode deeomposition definitions for tlie temperature perturbations, radiation 
dipole and the scalar Newtonian and curvature perturbations respectively we can write the 
mode decomposition of (148) and (149) as: 



R 



1 + R 



1 + Ra a 
kl 



ST « -H^A-^H + -^fi . 

a 3 



(152) 
(153) 



Upon ignoring the expansion coupled term (i.e the small scale approximation) and including 
the curvature fluctuations, since this makes the resulting equations applicable up to the Jeans 
length (above which the matter would not be gravitationally bound), and substituting the 
second equation (153) into (152) we obtain the well know equation describing the acoustic 
oscillations in the radiation [39]: 



6f" 



R' 

l + R 



ST' + k'^clST 



-<f>" 



R' 

1 + R 



H 



$4 

3 



(154) 



Here we have used conformal time (since we are now working in the conformal Newtonian 
frame) . 

6.2.5. The Dispersion Relations and Photon Damping Scale 

In this section we derive the dispersion relations for small scale anisotropies, where the focus 
is once again on developing generic covariant equations in parallel to the usual gauge invariant 
treatments [47, 51, 41]. To this end, we begin by iterating the baryon velocity equation in 
much the same manner as we iterated the integrated Boltzmann equations for the radiation. 

We begin with the the baryon relative velocity equation (131) which is once again inverted 
in order for it to take the form in (133). This is then turned into the basis of an iteration 
scheme in terms of the scattering time: 



V(n)i 



R 
k 



(155) 



Using this equation and the zero-th order tight-coupling approximation, i'"o)b ~ '^o, the co- 
variant second order baryon velocity equation is found: 



R 



i?2 



(H + H')v^ +Aa + HAa + 2Hv^ , (156) 



where we have dropped terms of 0{k~^). We now consider the small scale version of this equa- 
tion, by ignoring terms scaled by the Hubble parameter H and effects due to the gravitational 
potentials (see (66)): 



R r . . -, R"^ 

To - — [Ta + Aa\ 
AC 



2 ^'^ ■ 



(157) 



Expanding the transport equation for the second order radiation quadrupole (122) to first order 
and ignoring the shear contribution which is negligible in the almost- Friedmann-Lemaitre small 
scale limit, we obtain 



ab 



b) 



(158) 



2SWe use, as before, D^lnT « ^STQa , fa » fiQa , £>a*A = ^^AQa and T)a^H « ^^uQa- 
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where again we retain only first order terms. We now substitute (157) and (158) into the 
radiation dipole evolution equation (129) to find: 

-fa = {Aa + Da InT) - R{fa + Aa) + E?k-^fa - |k-1D^D<„t,) , (159) 

which can be compared to the mode equation A- 11 in Hu and Sugiyama [41]). This is the 
key equation from which we will now proceed to recover the dispersion relations and hence 
standard damping scale results. The key-point here is that the diffusion damping is second 
order in the scattering time, while the acoustic oscillations are first order. 

We now take a covariant mode expansion of the necessary quantities: Ta = TiQa, Aa = AkQa 
and Da InT = {^ST)Qa, together with the well known result (see Paper I [18] for the general 
case): D^Qab = —^{ak)~^{—k^ + SK)Qa. We consider only the flat case here, so we set K = 0. 
Putting these all into the second order radiation dipole equation (159) we obtain: 

-fi « (A + (f 5T)) + i?(fi +A) + R^k-'n + TE^-'^n . (160) 

Now we use the WKB approximation: 

Ti oc expi J {ijj/a)dt , (161) 
and drop terms scaled by H and a^a R. This gives 

+ R)^n « [(f <5T) + (1 + R)A] - R^k-'^n + ^f^-'^^n ■ (162) 



Now, in order to deal with the terms arising from Aa and In T we consider the covariant 
radiation monopole perturbation equation (15): 

(D„ InT) + ^BaO + H{Ba InT + Aa) ^ -^D„(Der^) . (163) 

We find, on taking a mode expansion and applying the WKB approximation again, dropping 
terms scaled by H and the expansion gradient (the small scale approximation described in 
section 2.5), that 

i-aST)^l^{-ico)n. (164) 

Upon substituting (164) into (162), factoring out the dipole coefficient ti, writing the differ- 
ential optical in terms of conformal time: k~^ = {aK')~^ and multiplying through by icoa we 
get: 

uj^{l + R) + R^K'~^uj^{iuj) - ^K'~^k^{iuj) « . (165) 
On re-arranging terms we finally obtain 

Splitting uj up into its natural frequency cjq and the diffusion damping term 7 and then solving 
the quadratic for the frequency w we obtain 

CO w ±coo + i'y , (167) 

where 
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Here kjj is the diffusion damping scale and c,, the baryotropic speed of sound in the matter. 
We get oscillations when k < ko (see the next section), and diffusion damping when k > kn, 
leading to a damping envelope. 

The above covariant result is equivalent to the analytic small scale approximation scheme 
of [39, 40, 41], who have demonstrated its robustness and closeness to more precise numerical 
studies^^. 

This derivation can be easily corrected to include the effect of anisotropic scattering (which 
breaks the tight-coupling approximation) and polarization. This is done by correcting the 
scattering terms in the calculation of the damping scale. Wc follow the approach of Kaiser 
[47, 39, 40, 41] and include the correction /2 to find the modified damping factor: 

« k^k^l , (169) 

where, first, for the anisotropic effect, f2 = jq, and second, to compensate for the polarization, 
/2 = |. Hence we get the diffusion damping envelope for n when k > ko- 

6.2.6. The Temperature Oscillations 

To examine the solution when k < kn, we take a mode expansion, using the Qa^'s defined 
in Paper I [18]: 

a{B^lnT)>' = k{ST)Qe ^ aDe(D»D„ InT) = (-:^)(^T)Qe , (170) 
where the driving term is written in terms of a generic potential defined by 

3 2 r,3 

-yDe(D2$^) « -a^{DSy + yDe(DM„) and a^B,{BHp) « +^^fQc ■ (171) 

From (145) we find the oscillator equation (once again working in the Newtonian frame cf. 
[40] eqn. (B3) and equation (154)). 

The sohition at first order must be convolved with the damping envelope, found from the 
dispersion relations, in order to include the damping cut-off. We have waited until the form 
of this damping is known, so that it can be easily included by simply replacing the natural 
oscillator frequency, uq, with w = coq + ^7, which now includes the damping. We use equation 
(145) in a manifestly gauge invariant form, to find upon mode expanding 

{6T)" + k^cl{6T)^-'^<^>F , (172) 

where the sound speed Cg is given by [3(1 + R)]~^- This gives the well known solution (see for 
example [39, 44] 

{ST){r]) « [{6T){0) + (1 + R)^f] cos(fcr«) - (1 + R)^f, (173) 

where the sound horizon scale is given by = / Csdrj « CgT]. This describes the source term 
for the acoustic oscillations with the isocurvature term dropped The last term on the right 

^^For numerical integrations of the covariant scalar temperature anisotropy equations we would refer the 
reader to [11]. 

^^The other part of this solution comes from sin(fers ) ; this is the isocurvature part ( ST' (0) ^ 0) , giving the 
full solution: 

{ST){r,) = [(fccs)-i(5r)'(0)] sinikrs) + [iST)iO} + (1 + R)^F]cos{krs) - (1 + R)^F. 
In this paper we consider adiabatic perturbations only so (i5T)'(0) = 0. 



JT+W 
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gives the Sachs- Wolfe effect, due to the potential; the other terms give the adiabatic acoustic 
oscillations. 

By putting this back into the momentum flux equation (136) in mode form we obtain 

t[ {v, fc) ^ - (kST) - ak^F , (174) 

and using the temperature anisotropy mode expansion that results when (126) holds, the 
Doppler contribution to the temperature anisotropics at decoupling can be found. These are 
given by 

Ti(r/, k) « (-A;-^)3[(alnT)(0) + (1 + i?)$„]c« sin(fcr«). (175) 



7. TEMPERATURE ANISOTROPIES FROM INTEGRAL SOLUTIONS 

In this section we derive the explicit form of the temperature anisotropics, using the general 
u"-frame integral solutions and the tight-coupling approximation solutions. 

Before wc do this, there are two important points that need to bo disciisscd. Firstly, where 
does the high-i? cut-off come from? The natural frequency of the system is loq and is set by 
the sound speed in the tight-coupled system to first order in Thompson scattering time - the 
oscillator equation. The diffusion damping is a second order effect whose correction is found 
by deriving the dispersion relations at second order in the Thompson expansion for the baryon 
and photon momentum equations and using the WKB approximation to find the new oscillator 
frequency to be : w « + il- The assumption that the anisotropics are sourced by these 
oscillations in the radiation is the key to the physics; the initial conditions before free-streaming 
are tai « for ^ > 2 (where the anisotropic correction at ^ = 2 is included as a correction to the 
damping scale). The free-streaming radiation transfer function (the spherical Bessel function 
in k-space) is then convolved with the initial power sourced by the oscillations in the average 
temperature and the dipole (essentially a cosine and sine function in k-space respectively). 
As free-streaming continues power is shifted up into the higher £'s, as the radiation transfer 
functions' maximum is near I oc kAr] (where Ar] is the elapsed time since last scattering). 
This maximum moves into higher £ for longer times to give one a sense for how the power is 
distributed through the i±l moment couplings in the almost-Friedmann-Lemaitre integrated 
Boltzmann equations. Obviously at high-^ (which corresponds to high-Zc) the amount of power 
surviving the (k/kn)^ damping will be very small, and hence as time progresses the peak in 
the transfer functions drops off for highcr-k, thus giving the high-^ cut-off. 

Given that in free-streaming there is no diffusion damping to cut-off the high-A; power, the 
truncation of the hierarchy is only consistent and meaningful if the significant anisotropy signal 
is sourced from the initial conditions at \ow-£ near to tight-coupling. 

Secondly, why do we find angular variations in the temperature anisotropics that are entirely 
due to oscillations in the dipole and monopole of the radiation - particularly given that in 
relativistic kinetic theory one has time- like integrations 

This is important for non-linear extensions of the almost-Friedmann-Lemaitre model [59], 
since care should be taken when treating ensemble averaged quantities integrated down the null- 
cone as being generally equivalent to those integrated down a single time-like world-line. The 
issue will become even more complex once the Gaussian averaging necessary for the construction 
of the angular correlation functions is relaxed. It is this, along with the weak Copernican 
principle that makes possible the extension of the data here and now, as integrated in a little 
box down a time-like world-line to last scattering, to global statements. 

^*Only in the almost- Priedmann-Lemaitre sub-case can one think of the averaged relationship between the 
null-projection 



1+3 COVARIANT COSMIC MICROWAVE BACKGROUND ANISOTROPIES II 



35 



Now we return to the issue of sourcing the temperature anisotropies from oscillations in the 
dipole and monopole. The integral solutions give the full almost-Friedmann-Lemaitre solution 
to the radiation anisotropies given the appropriate initial data. Generically it is the transfer of 
power from low-£ initial data up into high-f , with the accompanying reverse transfer of power; 
the £ zt 2, £ zL 1 transfer of power where there is a wall at £ = but none at high-^. There 
is strictly no £max unless the geometry is exactly Robertson- Walker [30]. The tight-coupling 
approximation gives the monopole and dipole at the ^ = initial wall. Diffusion damping in the 
initial power near last scattering as well as additional damping through slow decoupling cuts 
off the transfer of power. This power is sourced by those initial conditions and therefore give 
the temperature anisotropies in terms of the dipole and monopole alone. Additional integrated 
effects will only modify the primary projection, effectively leaving the cut-off unchanged. 

7.1. Temperature Anisotropies 

The integral solution gives the projection of these conditions at last scattering onto the 
current sky. These results arc covariant. 

The diffusion damping can be found from the dispersion relations arising from the coupled 
baryon-photon equations using the WKB method in tight-coupling which was described in 
section 6.2. We can now construct the primary source temperature anisotropies, by first con- 
sidering the acoustic and Doppler contributions in the free-streaming projection (94). Our 
solutions for Bq and Bi are first order in scattering time expansion hence to include the second 
order diffusion cut off, we multiply the solution through by exp7 (or exp7* if we want to 
include the corrections for the anisotropic and polarization effects). 

We can re- write 7 or 7* in terms of the damping scale, kn { cf [40] eqn (A7)-(A8)). 

7.1.1. Sources of Temperature Anisotropies 

In the manifestly gauge invariant integral solution of the temperature anisotropics in almost- 
Friedmann-Lemaitre spacetime, using the slow decoupling approximation with scalar perturba- 
tions (98,102), there is additional complexity of having to deal with shear terms. We therefore 
choose the Newtonian frame in which the shear contribution does not appear^^. We are then 
able to recover the canonical scalar treatment [81, 75, 35, 39, 40]. This can be written out in 
terms of the scalar perturbation potentials and the temperature perturbation, ST, the 
radiation dipole n (through slow decoupling) and the baryon relative velocity, Vb, all within 
the covariant approach (including Doppler contributions (111)): 



fo^l^^il ^Sp{k,v*).u(kAv*)+ [ dr]{SDisw + Sisw)je{kAT]) . 

+ IJ ' « ' Jri, 



(176) 



primary sources 



secondary sources 



Here the source terms are given by: 



Sp{7j,k) = V{r,o,k) {6T + <Pa)+h'vb 



-kfi + {n'v'g + k"vb) 



SDisw{ri,k) = Ve-^^/'^^) 
Sisw{v,k) = Ve-('=/'=^)' [^'A-^'H]{k,v)-aH[6T + 3^A]{k,v) 



(177) 
(178) 
(179) 



In the first line, on the RHS, we have the Sachs- Wolfe and acoustic sourced projection effects 
which are the primary sources. The next term describe secondary Doppler effects during slow 



^^Incidently this would also remove any problems we may have with the introduction of high-^ truncation as 
discussed in Appendix E 
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decoupling. The last term models the integrated Sachs- Wolfe (ISW), the late-ISW and early- 
ISW effects respectively. We have used Atj* = rjo — rj* (the conformal time difference between 
the surface of last scattering and the time of reception, here and now) and Ary = ijq — rj. 

In this way we have recovered the solution from the exact equations in a systematic manner 
using the 1+3 CGI formalism rather than recovering the solution from a perturbation theory 
about a foliation of Robertson- Walker surfaces of homogeneity. Note that there is no monopole 
temperature anisotropy in the CGI approach while there is one in the canonical treatment. 

The angular correlation functions, C^, for the general small temperature anisotropy case have 
been derived in terms of a superposition of homogeneous and isotropic Gaussian random fields 
with respect to the temperature anisotropy multipoles to find the multipole mean-squares, 
(ta<,t'^*), in terms of the an ensemble average (see Paper I [18]). The relationship between 
the multipole mean-squares (in the almost- Fricdmann-Lcmaitrc case) and the mode coefficient 
mean-squares |r^p in terms of the Robertson- Walker mode functions Ge[Q] are also given in 
Paper I [18]. We do not discuss these further here. It must be emphasized that the multi- 
pole mean-squares are given for general geometries, while the mode means sqiiares are only 
for almost-Robcrtson- Walker geometries. This is why the application of the non-linear exten- 
sion is possible; where the corrections to the almost-Friedmann-Lemaitre standard model are 
calculated using the multipole formalism of [59]. 

7.1.2. Sachs-Wolfe Effect and the Acoustic Source 

In the Newtonian frame (177) we find that the Sachs- Wolfe effect arises from a combination 
of the DalnT + f (^^ + ^a)<3o- From the solution to the oscillator equation in the 

Newtonian frame (eqn. (154)) along with ^ CsTj and ignoring the time evolution of the 
potential (for the equivalent canonical version see [39]) we find 

Sf{r], k) + $a(??*, fc) « [5f (0, k) + {l + i?)$A(0, k)] cos(fcrs) - i?$A(?7*, k). (180) 

Upon taking the matter dominated limit {R « 0(e^)) and then using the adiabatic assumption 

in tight-coupling (5T(0, fc) \/^{Q,k) ^ -|$a(0,A;) (adiabatic flat CDM model) and finally 
taking r* ~ (by using that fcry* <C 0) we recover the usual results as in [39, 43] 

[6f + ^a]{v*, k) ~ +1^a{0, k){l + 3i?*) cos(A;r:) + R,^a{0, k) ~ 1^a{0, k). (181) 

Here we have used that ^^(j?*,^) ~ ^"^(0, fc) from the Einstein de-Sitter result that the 
potential is constant if we drop the decaying mode. The physics of the Sachs- Wolfe effect as 
opposed to the acoustic oscillations is then quite clear. There are no oscillations just an imprint 
due to an acceleration potential. Although the adiabatic assumption is invariant to order 0{e), 
the relationship between the Electric part of the Weyl tensor and the perturbations are not - 
so we will always use the matter perturbation in the total frame where it can be easily related 
to the Newtonian potential and the temperature perturbation in the Newtonian frame where 
the oscillator equation takes on a useful form. In a similar manner we find the explicit form of 
the radiation dipole and acoustic oscillations in the slow decoupling era: 

5f{k, 7j) « i$A(0, k){l + 3R) cos{kcsv) + (1 + R)<^a{0, k), 

fi{k, 7?) « $a(0, fc)(l + 3R)cs sin(fcc,7?). (182) 

7.1.3. Weak-coupling for the Small Scale Solution 

Here we briefly consider the integrated Sachs- Wolfe effect in connection with the weak- 
coupling approximation. The idea is that the anisotropics fall with i more rapidly than a 
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simple projection would imply. What one has in mind is the situation where the anisotropy 
contributes across many wavelengths of the fluctuation allowing cancellations on small scales; 
the secondary sources, in particular the term — ^'h), varies slowly on small scales [43]. 
Specifically we consider the situation where 

dr?Ve-('=/*'-)^ |($^ _ $^)} j,(fcAr?) « ,J ~{^'a ' ^'hK^V^MKvo) , (183) 

where wc used f^^ je{kAr])dr] = [^/2k][T{{e+l)/2)/r{{£+2)/2)]. The weak-coupling solution 
is implied by the assumption that: 

{<i>H - ^a) < H^H - ^a) . (184) 

This is nothing more than a useful approximation allowing the direct construction of analytic 
solutions. Some care should be taken when using the weak-coupling approximation when trying 
to estimate the early-ISW (the free-streaming analogue of the acoustic driving effect) and 
late-ISW (due to non-vanishing curvature or cosmological constant, which will dominate the 
expansion rate at latter time) effects. Ideally one should evaluate the slowly-varying function, 
which has been taken out of the integral, at the i-th peak; r]e = r]o — {i+ \)/k rather than at 
Ary* = % - V*- 

In the case of the early-ISW effect, since it satisfies neither the tight-coupling nor weak- 
coupling criteria, our approximations schemes here are not entirely appropriate; its decay time 
and wavelength are comparable [39, 40, 43]. However we can use the weak-coupling in the 
case of the late-ISW effect because cancellation effects lead to damping on small scales. The 
temperature perturbation decays, and hence the potential decays on the order of the expansion 
time near the end of the matter dominated era. The photons will free-stream across many 
wavelengths of the perturbations below the Hubble scale, leading to cancellation and damping 
effects. 

The important point about the late-ISW effect is that there will be an imprint due to the 
exit from the matter dominated era into a A or curvature dominated one. In the context of the 
equations here, we can consider A-dominated effects by investigating the evolution equations 
for the potentials in a A-dominated case. 

In this case the effective expansion changes, given here for late-times (well after matter- 
radiation equality): 

H'^ = a-^floH^ + \a (185) 

for the r^o + ^^A = +1 model. This has however been dealt with in depth in the literature (see 
for example [43, 40]). 

7.1.4. The Mode Coefficients 

The angular correlation functions, measured here and now (xq) can then be computed using 
the the mode coefficients from Paper I [18]: 

The final step is to construct the temperature anisotropy solutions in terms of the matter 
power spectra. We discuss this in the next section. At linear order, there is no important 
difference in our solutions from those found in the canonical treatment [39, 40], however the 
advantage is that we can easily relate our formulation of the temperature anisotropics to the 



(186) 



38 



GEBBIE, DUNSBY AND ELLIS 



mean-squares of the multipoles (ta^ ) and hcnco to the mean-square of generahzed temperature 
anisotropies, (H^^) [59]. This was not attainable in the canonical treatment. 

7.1.5. The Multipole Coefficients 

The mean-square of the multipole moments are related to the almost-Friedmann-Lemaitre 
mode coefficients by (see Paper I [18]) 

{ta,t^') « j k^dk\n{k,r]o)\^ , (187) 

giving the angular correlation function (see Paper I [18]): 

d = Ai{2e+1)-' {ta,t^') . (188) 

This relates the matter fluctuation amplitude at last scattering to the present day via the 
Newtonian like potential, which in turn may be related to the matter power spectrum directly. 
A plethora of numerical studies of Doppler peak features exist in the gauge-invariant literature 
(see for example [42] and [74]). 

We will now summarize the standard picture of acoustic peak characteristics before discussing 
the matter power spectrum. 

7.2. Some 'Acoustic peak' Characteristics 

The key features of the standard model of cosmic background radiation primary sourced 
Doppler peaks are listed below [39, 40, 43, 44, 79, 46]. The standard model of acoustic peak 
formation has been given in its analytic form above. 

1. The j — th peak positions (as given in the flat adiabatic case) is given by 



£j w kj \re{v)\ = 



jn , (189) 



where re is the comoving angular diameter and is the sound-horizon near decoupling. Notice 
that je{kAr]) peaks near i ~ kArj, where free-streaming projects this physical scale onto the 
angular scale 9Arj on the current sky. The first peak is dependent on the driving force which 
is model independent. It provides a way of fixing the angular diameter distance - when using 
the almost-Friedmann-Lemaitre assumptions. 

2. The acoustic relative peak spacings are given by 

kj+i -kj=kA^— Ai^eA = kAre ■ (190) 

The peak spacing is fixed by the natural frequency of the oscillator: w = kCg, which is inde- 
pendent of the driving force. The factor is the photon-baryon sound speed. 

3. The peak ratios arise from the angular power spectrum ratio Cjk = Cj/Ck found from 
in terms of jr^p or {ta^t^'^)- 

4. The damping tail provides yet another angular diameter distance test of curvature [43] via 
the damping scale kjj . The peak spacings A(. and the damping tail location Id depend only on 
the background quantities, they are robust to model changes, assuming that secondary effects 
do not overwhelm the signal. The diffusion scale near decoupling is the angular scale of the 
wavenumber ku ^ ^J k/Arj --^ y/ajArjh ~ Ar]~^. Given that the damping tail of the acoustic 
oscillations takes on the form e"*^*^/*^"^ , ko can be found from (as derived previously): 
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for no anisotropic stress [43, 47]. This can be used to find that the damping tail location is 
(■D = kore. 

5. The damping tail position to peak scale, £d/^a, is a good measure of the number of peaks 
and gives an indication of the delay in recombination independently of the area distance. 



8. THE MATTER POWER SPECTRUM 

In this section wc related the temperature anisotropics in the case of the Sachs- Wolfe effect 
to the matter power spectrum using the 1+3 CGI formalism. We also deal with basics of the 
normalization of the cosmic background radiation angular power spectrum to the matter power 
spectrum. 

Up to this point everything we have derived has been found covariantly from general relativity 
and from relativistic kinetic theory, however on small scales we do not have a covariant analytic 
derivation of the transfer function, although it is well known for large scales, so this is the only 

gap in our treatment 

The aim here is to try to predict the matter distribution today from the cosmic background 
radiation spectrum. Prom the cosmic background radiation spectrum one finds the matter 
power spectrum (as fixed by the cosmic background radiation data today) from which the 
current distribution of matter is found (in terms of as). 

Briefly, in the literature there have been two methods for normalizing the cosmic back- 
ground radiation data in relation to COBE: firstly the a{10^) normalization where the r.m.s 
temperature fluctuation was averaged over a 10^ FWHM beam, and secondly the Qrms-ps^^ 
normalization which uses the best fitting amplitude for a n = 1 Harrison-Zel'dovich (HZ) spec- 
trum quoted for the quadrupole, (Q). The first method had the advantage of being a model 
independent way of fitting the data - it is observationally determined, however, because of 
this it does not provide the most accurate normalization for a specific model and care must be 
taken in order to properly deal with cosmic variance. The second approach is model dependent 
and works well only for the HZ spectrum [20, 9]. 

Improved, more general normalization schemes have now been developed and there currently 
seem to be two favoured ones, both having variations of the HZ spectrum in mind. It is these 
two schemes we now briefly consider here. We emphasis the so called CDM-models with 
CIq = 1. These models have a baryon fraction ^Ib with the rest of the matter being made 
up of massive Cold Dark Matter (CDM). Of these the most easily dealt with is the so called 
standard- CDM model, where initial fluctuations are assumed to have a Gaussian distribution 
and have adiabatic, scalar density fluctuations with a HZ spectrum on large scales: 

P{k) a fc"-^ , n = 1 , Ho = SOkm-^Mpc"^ , Os = 0.05 . (192) 

There are three other popular models: the Adiabatic CDM model of Peebles, which is a 

version of standard-CDM; the /CDM model (which is an isocurvature version of the Adiabatic 
CDM model; and lastly the ACMD model which is the large cosmological constant version of 
standard-CDM. More recently Hu has introduce a Generalized Dark Matter model (GDM) [37]. 

^'^ There is still some hope that the almost Priedmann-Lemaltre Quasi- Newtonian models, will be useful in 
this regard [31] providing us a covariant derivation of the Harrison-Zel'dovich spectrum. 
^"^Qrms-PS « ^Ca and Cf"- ~ 
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8.1. The Power Spectrum 

The Power Spectrum V{k) is fully specified by a shape and normalization (see Appendix 
J for additional details and references). The normalization is fixed by the amplitude of the 
temperature fluctuations in the cosmic background radiation^^. 

There are two normalization scales, large and small (see Appendix J). On small scales the 
normalization is expressed in terms of cis, the ratio of the r.m.s mass fluctuations to the galaxy 
number fluctuations - both averaged over randomly located spheres of radius 8h~^Mpc; this is 
the variance of the density field in these spheres. 

On large scales the power spectrum is found by fixing the shape function [19] which is a 
measure of the horizon scale at matter radiation equality. We look at these now for the case 
of a fiat matter dominated almost-Priedmann-Lemaitre model (the open case is considered in 
Appendix K). 

In the 1+3 CGI approach the emphasis is on the quantity aDalnp^^^ « aS^^jQa- Here S^^^ = 
aS^{0,k) and {d^{k,0)S^(k' ,0)) = {2T:fV{k)5{k-k')/k^ for |(5„(fc,0)|2 = V(k). It becomes 
preferable in the context here to use the variable A: ^A(fc,t) = (5„(fc,t)^^ . We will, however, 
be using the power spectrum for the gauge invariant acceleration potential because we consider 
the situation of matter dominated adiabatic perturbations. In the general situation it is best 
to use Note, in the matter dominated situation we can use that = —^h and write 
everything in terms of either or ^a- 

The relationship between the power spectrum of the acceleration potential and the power 
spectrum of the matter perturbation is 

\^Aik,0)\^ = P^^{k), and |A(fc,0)|2 = P(A:), (193) 

where these are then related to each other in the total frame via the constraints for the electric 
part of the Weyl tensor for the fiat case: 

P^Ak,0) = ilH^no^f k-''P{k). (194) 

— .2 

Here $ = —^^a and as noted before: 

i^pZ'){k,t) « -{3H^no)-'D{t)[kHA{k,0)]. (195) 
and D K a K for a fiat almost-Fricdmann-Lcmaitre dust model. We will also use ao = +1. 

8.2. Relating the Power Spectrum to Temperature Anisotropies 

We now combine the CGI Sachs- Wolfe effect (due to potential fiuctuations) and the mat- 
ter power spectrum to express the anisotropies in terms of the matter power spectrum near 
decoupling. 

8.2.1. Flat Matter Dominated Almost-Friedmann-Lemaitre 

We first consider the fiat almost-Friedmann-Lemaitre model. Using (193) and (181, 177 and 
176) (that is with K = 0) we obtain 

^1^^ ~ +^^A{k,vMkArj) , (196) 

Here we find the angular power spectrum from (186), (196) and (194): 

Cf- = ^ {H^^o^y \ J §P(fc).?1(fcAr?*) . (197) 

^■^It seems that models normalized to COBE which are fixed in both the amphtude at small scales and by 
the shape predicted by standard- CM D, are inconsistent with observations [49, 9]. 
33A(fe)2«.f^|5„(fe)p = |^Afc»+i 
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where P{k) is given by either (J.l) or (J. 3). Now if use that (Oo-D*/a*)^ ~ f^o'^^ we find the 
standard result [20]: 

1 f rlh 

Cf- = ^ffo'^o j ■j^P{k)jj{k^rj.) . (198) 

The next step is to normalize the matter power spectrum to any given structure formation 
theory on both small and large scales, this is done in the next section in the case of the 
standard- CDM model. 

8.3. Setting the CDM Normalization 

Prom the previous section we can now relate the angular power spectrum for the matter 
dominated flat (K = 0) almost-Friedmann-Lemaitre model to the current matter power spec- 
trum on small and large scales using the two different normalizations to the standard-CDM 
model. In what follows we will use the useful result: 

Jo z^^' ~ 2^+2 (m/2)! (^ + f + 1 )! • ^'^^^ 



8.3.1. Large Scales 
Prom (198) and P{k) = Ak''-'^ (J.l) we obtain 



C, = ^AH^nl-'' J ^k-'jf{kx) , (200) 

and on using (199) for m = 2 (n = 1) we find that 

A 1 
C, = 2^fo^^o''(2£^3)(2£+l)(2^-l) ■ 

Of more immediate interest is the angular correlation function for matter below the horizon 
scale near decoupling, as this is what is used to normalize the angular correlation function. 

8.3.2. Small Scales 

Por small scales one can use the erg normalization via (198) and P{k) = BT'^{k)k"' (J. 3) and 
the parametrized transfer function (J. 5). This gives 

C, ^ ^H^i^l-^' J ^Bk^T\k)jUkAri.) . (202) 

For example using the standard-CDM m.odel with n = 1, h = 0.5, fls = 0.05 and F = 0.48 we 
can relate B and A to the scale of fluctuations near horizon crossing and the quadrupolc: 

B = 2Tr^A = (677^/5 )(g)/r2 . (203) 

To normalize to as we can choose as — 1.3. This then allows us, in principle, to invert the 
angular correlation function to find the matter power spectrum once the initial spectrum is 
known. 

We are however more concerned with normalizing the radiation angular correlation function 

to CDM on horizon scales near decoupling. In this regard we once again begin with the matter 
power spectrum (we set T{k) = 1) and on using (199) for m = 1 we find that: 



^{CDM) ^ 



-H,n, B— 



^ (204) 
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8.3.3. Ci Normalized to Adiabatic CDM 

Finally we normalize the angular correlation function Ci (found from 102) to the potential 
fluctuations normalized to Adiabatic CDM (204) above: 

-1 

£{£ + l)Ce . (205) 

The angular correlation function Ci can be found from the primary (177) and secondary (the 

( p) ( s) 

integrated part of 102) source that make up the total: C^ = C\ -\- C\ . This then allows 
one to remove that part of the angular correlation function arising from the standard Sachs- 
Wolfe potential fluctuations leaving a signal which is dominated by the photon primary and 
secondary sourced physics. The convention is to use rather than Q [46], so we have from 
(198 and 204) for the flat case: 

De = t{i + 1)Q = y |^p*(fc)|i,(fcx)P using V*{k) = j , (206) 

where kg is the normalization scale, a gives the deviation from the power law, and rig give the 
scalar power law index. The angular power per ln£ is + l)/47r)C^. 



2it' 



9. CONCLUSIONS 

In this paper we have carried out a covariant analytic time-like integration reproducing 
the well known primary effects generating the 'acoustic' peaks measured here and now for 
an almost-Friedmann-Lcmaitrc universe with adiabatic sc;alar perturbations. We have also 
demonstrated how, in the CGI formalism, the angular correlation functions are constructed in 
terms of the matter power spectrum and normalized on large and small scales for standard- 
CDM, given appropriate approximations for the transfer functions. 

As pointed out initially, the aim of this paper was to clarify the link between the standard 
gauge- invariant and CGI treatments of cosmic background radiation anisotropics and provide a 
strong basis from which to tackle non-linear and gravitational wave effects using CGI methods. 

Some of the key outstanding issues are: 

1. How does one deal with anisotropic scattering and anisotropic stresses before and during 
decoupling within the covariant approach, speciflcally in such a manner so that consistency 
is maintained when using general relativity, its covariant linearisations and relativistic kinetic 
theory. 

2. The small anisotropy equations developed in [59] with the application to spacetimes 
with arbitrary anisotropy and inhomogeneity have yet to be properly investigated; these be- 
come applicable when one ignores the Copcrnican principle that underlies the almost-Ehlers- 
Geren-Sachs theorem, which in turn provides the theoretical basis for using almost-Friedmann- 
Lemaitre spacetime dynamics. An investigation of their consequences on the cosmic background 
radiation may provide an alternative method of testing the Copernican principle other than 
the Sunyeav-Zel'dovich effect or via the use of polarization maps. 

3. There is a need to find a working non-Gaussian treatment from which one can construct a 
generic characterization of observables here and now (other than the angular power spectrum 
alone (see Paper I [18]) and, second, finding an ab initio covariant analysis of transfer functions 
extending the post-Newtonian treatments which use periodic boundary conditions. 
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In the next paper in the series, Paper III [15], we hope to establish, in a complete fashion, the 
relationship between the null-cone integrations (favoured in the literature) and the time-like 
integrations (found in the exact relativistic kinetic theory). 
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APPENDIX A 
Some Useful Almost-Friedmann-Lemaitre relations 
Some useful covariant linearised differential identities are given below [60]: 

curlD„V = -2V'Wa , (A.l) 

(D„V) « D„ (dV) + i (p - 3//') D„V + 2V'curla;„ , (A.2) 

(D„V)' « D^V- - HBa^ + Ma , (A.3) 

(aD^J^-A,)' ~ aD^j^-^, , (A.4) 

(D^T/;)' w D^?/; - 2i7DV + V-D^^a , (A.5) 

D[aD(,]yc « I (31/2 _ ^) y^^,^^^^ ^ (A.6) 

D[„D(,]5^'' « I {3H^ - p) 5[„('=/i6]'') , (A.7) 

D^curlV; « 0, (A.8) 

D'-curlS^b « icurl (o'-^.b) , (A.9) 

curl curl « D„ (d^h) - DV„ + § (p - 3if2) , (A.IO) 

curlcurl5„6 « |D(„D'=5b)e - D^Sab + - SF^) , (A.ll) 

£)''D''D(„H> « \'D'^dwV +{p-iH'^)dwV , (A.12) 
where the vectors and tensors are 0(e) and Sab = S'(ab). 

APPENDIX B 
Integrated Boltzmann Equation Relations 

Here we repeat some useful results from [59]. The integrated Boltzmann equation is : 

/•oc POD r r) f "1 

L{f)E^dE = J^ E^dE^p''daf-r\yp^-^^=j^ E^dE. (B.l) 

Here f is the single particle distribution function, C[f] the scattering correction. The Liouvile 
operator for the photons, L{f) = ^ and d/dv = p°'Va is the null derivative. Hence we find 

L(/(i.,.Ne")) = |=p"V./ + ^|^, (B.2) 

as in [59], so that given the identity (see [59]) 

^ = -E"" [\@ + A„e« + Gabe^e^] , and p» = E{u'' + e") , (B.3) 

we have covariant derivation of the almost-Priedmann-Lemaitre integrated Boltzmann equa- 
tions. 

APPENDIX C 

Scattering Strength Expansion 

The almost-Friedmann-Lemaitre integrated Boltzmann equation for Thompson scattering is 

B + f + e'^DaT ~ t-^ [v^ea - r] . (C.l) 
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This enables us to find 

t{x\ e") = v'^Ca -tc[B + f + e«D„T] . 
We now systematically approximate (C.2) in terms of the smallness parameter tc 

T{n) « V'ea -tc[B + f{n - 1) + e"Dar(„_i)] . 
Up to second order the following anisotropies are recovered: 
T(o)(a;',e") w w^e^ , 



B- 



(C.2) 

(C.3) 

(C.4) 
(C.5) 

(C.6) 



In order to carry out the solid angle integration over the sky (116), the following results 
will be useful. From the normalization of the direction vectors, e", along with the recursive 
definition of O"^' , it can be demonstrated that 



s'^e^e^DaDbVc = O^^^DaDbVc - ^ ( 0"D(„D(,)u'' + O^DaD"vb) . 

o 



Using the orthogonality conditions we find, first, from (68): 



/ 



I4n 

and second, 



47r 








+ si 




.3 







(5^ 



15' 



/ 



(C.7) 
(C.8) 
(C.9) 

(C.IO) 
(C.ll) 



The integration over the solid angle can now carried out resulting in an equation for the gradient 
of the radiation flux. 



APPENDIX D 
Integral Solutions 

We use the following notation: Ar^* = rjo — r]^:, Arj = rjo — rj and 6r] = r] — rj' such that 
5r]* =rj- rjt, and 5r]o =r] — Vo- 

Beginning with an integral ansatz of the form: 



Jo 



C,{n')Tf\5v) + Ci(V)^rr(M + C2(r;')-^^r('5^) 



92 



drj' 



(D.l) 



Using the Leibnitz rule for differentiation of integrals we obtain: 



drl 



Q or] I 
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+ 



Coiv)r^'\o) + Ci(,?)rf' (0) + C2{v)4'^" {0) 



If we hold rj' constant in the partial derivatives, it follows from (86) and (D.l) that 



(2£+l)(2£ + 3) 



d_ 



This gives us 
k 



(2^+l)(2^ + 3) 



Putting (D.l), (D.2) and (D.4) together, we find : 



rf + k 



(2£+l)(2^ + 3) 



- rg_i 



Co{r,)Tf\Q) + Cr{r,)Tr (0) + C2(7?)Tr (0) .(D.5) 



JO)' 



(0)' 



(D.2) 



(D.3) 



(D.4) 



APPENDIX E 
Truncation Conditions 

The Ellis- Trcciokas-Matravcrs treatment [30] makes various general statements about exact 
relativistic kinetic theory in the free streaming case. Of these probably the most important 
are: 

Firstly, if any four consecutive harmonics vanish, say those with I = Z/ + l,i + 2,Z/ + 3, i + 4, 
but those for / = L are non-zero, then 



= 



(E.l) 



This means, as Fai is non-zero, that the shear must vanish exactly: Cab = 0. This results 
arises from the requirement that lim^j^oo ^A, = 0. 

Secondly, it can also be shown that if the first 3 multipole harmonics are zero, i.e. , 
I = 1,2,3, once again the shear is necessarily zero: aab = 0: 



fab 



(E.2) 



Thus in both cases the resulting spacetimes are highly restricted, and do not include generic 
perturbations. 

A simplistic approach to linearisation will suggest that these relations and their implications 
can be ignored in (linearised) almost-Friedmann-Lemaitre universes, because the equations 
leading to these results are second-order relations and so can be dropped when linearising. 
However that argument is not correct, if one carries out a careful linearisation procedure: 
indeed both these statements will hold in almost-Priedmann-Lemaitre universes also. This can 
be seen as follows: although both Fai and (Jab are at most 0[1] (or 0(e)) in the almost-Ehlers- 
Geren-Sachs sense [73] , there are no zero or first order terms in the relevant equations leading 
to the above results, to explicitly linearise with respect to. Thus they cannot be dropped 
relative to larger (first order) terms in these equations, as there are no such terms; the first 
non-zero terms are second-order, and hence these equations with these terms must be obeyed 
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even if wc carry out a (first-order) linearisation. Thus they are both at most 0[2] equations, 
but are both still valid in the almost- Friedmann-Lemaitre universes^. 

What this means is that one must be very careful about any kind of truncation in the 
multipole hierarchy, even in the almost-Priedmann-Lemaitre universes - this includes not only 
the free-streaming case but the case with a generalized Krook equation for the scattering term 
[30]. What should be of interest is that in the matter dominated almost-Friedmann-Lemaitre 
models with scalar perturbations, any truncation leading to zero shear would suppress the 
perturbations, reducing the dynamics to that of an exact Friedmann-Lemaitre model. This 
doesn't mean that one cannot consistently damp the higher moments out, it just means that 
they cannot be formally truncated - when higher moments are ignored the consistency condition 
in the exact theory should still be taken into account. Where this makes a difference is, for 
example, on intermediate scales where one is tempted to drop everything with Z > 3, i.e. , 
when one is close to tight-coupling. However, this is, strictly speaking a truncation and hence 
problematic. This is why a perturbative analysis in the Thompson scattering time is important; 
one can in this way consistently, without truncation, build up the entire multipole divergence 
equation hierarchy perturbatively - provided one has a meaningful sense of smallness, in this 
case the relaxation time. We do this in the case of small scales, in this way decoupling a subset 
of the hierarchy from the full set of multipole divergence equations, as in the gauge- invariant 
formulation of Hu-Sugiyama [39, 40]. 



APPENDIX F 
Linking Different Expansions 
The gauge invariant and covariant mode expansion (6) and (7) [18, 11] in the almost-flat- 
Friedmann-Lemaitre case {K = 0) give the mode coefficient recursion relations for i>2: 



k 

-n ~ - 
a 



.(2f + 3)(2£+l) 
Multiplying through by fii, one finds 



(^ + 1) 
{2£ + 3) 



{21 - I) 



(F.l) 



(F.2) 



Changing from the proper time derivative in (F.2) to the conformal time derivative ' using 
dt = adrj we find 



-iPm)' ~ k 



(2£ + 3) 



(2£-l) 



(F.3) 



This can be immediately seen to be the same mode equation for £ > 2 as in Hu & Sugiyama 
[39] and Wilson & Silk [81] (eqn 7). 

If now I3i is replaced with aj^{2l + 1) = we find that (F.3) can be rewritten (on first 
multiplying through by (2^ + 1)~^) as 



-{2t+l){a-^n)' « fc [(£ + l)(a7+\r,+i) - ^(a7_\r,_i)] 



(F.4) 



This can be recognized as the form of the (. > 2 free-streaming integrated Boltzmann equation 
of Ma & Bertschinger [56] (c/. mode equation, 49 or 50) or Seljak & Zaldarriaga [69] 

equation 3d. Once again, it may be useful to remind the reader of our nomenclature, are 



^Although O(e^) <C 0{e), this doesn't mean that O(e^) = on its own, even though formally the notation 
O(e^) ~ is often adopted. These are not equivalent. 
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the mode coefficients (from tlie mode expansion) wfiile ta, arc tlie multipole coefficients (from 
tfie multipole expansion). This distinction is not made in the Bardeen- variable gauge-invariant 
treatments based on the Friedmann-Lemaitre mode expansion. 

The solution to the mode coefficients with respect to the time-like integration in the flat 
{K = 0) almost-Fricdmann-Lcmaitrc case, arc spherical Bcssel functions. This should not be 
surprising given that the recursion relation in the linear Friedmann-Lemaitre case is merely a 
projection from an initial section onto a sphere around here and now, modified as a result of 
Robertson- Walker expansion. 

APPENDIX G 
Almost-Ftiedmann-LemaJtre Equations 

Here some of the results arising from the Einstein field equations are listed, these are pre- 
sented in [59]. 

G.l. MAIN ALMOST-FRIEDMANN-LEMAiTRE EQUATIONS 

Wc give the almost- Friedmann-Lemaitre constraint, propagation and perturbation equations 
[57, 28]. Using the 1 + 3 covariant dynamical equations [32] [24], in geomctrizcd units : c = 
8nG = 1. With the projection tensors U% = -u'^Ub, h\ = 5\ + u^Wb, /i^'^''^ = and = 
DaC^bc = Da/i(,c, where the totally projected spatial derivative is D^,. The covariant derivative 
of is 

VfcMa = -UaAb + DaMb := -UaAb + (Jah + Hhab + ^abc'^" , (G-1) 

where we have written the vorticity in terms of the vorticity vector, Sabc = 'Habcdu'^ , and the 
shear Gab = <7<ab> ■ We consider a cosmological model with matter (dust) and radiation. Then 
the stress-tensor given by 

= (pM + Pij)w"w^ + 25("w'') + w"'' + ph''^ , (G.2) 

where energy densities pM and pn are the energy densities of the matter and radiation respec- 
tively, q"^ is the total energy flux, TTab is the anisotropic pressure of the radiation. 

We linearise about a Robertson- Walker model as explained previously (see section 3.2 and 
[25]). The convention A = B implies that A equals B in the exact theory, while A :^ B 
indicates that A equals B to at least 0[1] in the almost-Friedmann-Lemaitre sense, and lastly 
A Pi B \s retained to indicate that A equals B in relation to some other specific smallness 
parameter such as the relaxation time or ratio of radiation density to matter density. The 
almost-Friedmann-Lemaitre Gauss-Codacci relation (Hamiltonian constraint) is 

3i?~2/>+^e2. (G.3) 

Using the Ricci identities, the linearised propagation equations are 

Q + \e^-{diYA) + \{p + ip) ^ -\b, (G.4) 

W(a)+'2'Hw(^a) + \cm\Aa ~ 0, (G.5) 

(^{ab) + "^Ha^ab) - D(o^6) + -E^(ab) — 1^''^ {ab) ) (G.6) 
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and the constraint equations are given by 

(diva)" - ^D"e - curlw" ~ q'' , (G.7) 
Daw" ~ 0, (G.8) 
H {ab)- curl a (^ab)-'D (a i^b) 0. (G.9) 

Using the second Bianchi identities, the hnearised propagation equations are 

p+{p + p)Q ~ -SFB-DV, (G.IO) 
(p + p)A^ + DaP -D„S - 9„ - 4Hq^ - B^^b , (G.ll) 



^{ab} ^ 3^^(a6> _ curli?<°''> + ^ (p + ^ "^'^^"''^ ~ ^ff7r<"''> - iD<"g''> ,(G.12) 

ijM) + 3ijij(«6> ^ ~ ^curl7r<''''> , (G.13) 
and the constraint equations are: 

(divS), - ^D.p ^ - ^B\ab , (G.14) 

(div (p + p)a;a ~ -^curl^a . (G.15) 
The background Friedmann-Lemaitre equations are: 

p = 3H^ + ^, (G.16) 

p = -3if(p + p) . (G.18) 



G.2. SOURCE TERMS 

Prom [59], the source terms in the almost-Priedmann-Lemaitre multipole divergence equa- 
tions are: 

Bo{x) ^ -^D^Ta{x) , Ba{x)^DalnT{x)+Aa{x), Bab{x) aab{x) . (G.19) 

In the case of scalar perturbations in a matter dominated almost- Friedmann-Lemaitre universe 
[59], using the mode expansion [18], the electric part of the Wcyl tensor Eat (gravitational tidal 
effects), the acceleration Aa and anisotropic pressure nab can be related to their corresponding 
potentials: 



Eab = Hk,t)Qab, Aa = ^u{k,t)Qa, TTab = K{k,t)Qab ■ (G.20) 
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Using the Mode functions Q, Qa and Q(^ab), taking care of the additional wavenumber factors 
which get introduced ^ [18], the following decomposition will be used: 



(Dalnp„)*^ = S'^'{k,t)Qa = +^A{k,t)Qa , 



It then follows that ^ 



Boik,t) ^ ~-n{k,t) 

O Qi 



B2{k,t) ^ -2\p-\t)^k,t)\ 



(G.21) 
(G.22) 
(G.23) 

(G.24) 

(G.25) 
(G.26) 



In the case of matter domination we then find the relationship between the Newtonian potential 
and the matter density gradients: 



(G.27) 



Replacing a^p^ in terms of the density parameter and the curvature constant from the Freid- 
mann equation (See Appendix A) for both the K = Q (Oq = +1) and K <Q (Q.q < 1) matter 
dominated cases respectively, we obtain 



$(fc,i) « -(ifo'f^o)^(?7)A(fc,r?o) 



Hk,t) 



3k 



{k^ - 3K) (fio - 1) 

1 Kflo 



2 a (A;2 - 3K) (fio - 1) 



D{r,)A{k,Tjo) 



(G.28) 



(G.29) 



These are then the equations that will be used for the slow decoupling and free-streaming 
eras. 



G.3. THE TEMPERATURE MONOPOLE 

Using the energy conservation equations {£ = multipole divergence equation) ^ 

(D„ InT) + H{Ba InT + A„) + ^D^O ~ +iD„(DVe) , 



(G.30) 



^Once again one should not confuse the Fourier coefficients of the potentials here with those defined from 

Eab = D(a6>*(a;0' = Da$„(x*) and IVab = G{ab)^n{x*). 

^To see where these factors come from notice that a^'D'^DacQ = {—k"^ + 2K)DaQ =^ a'^D'^D/^^^^Q = §(— fc^ + 
3K)DaQ (after removing the trace) for the general I — th order relations see [18]. {—^)^^^(Ai)Q = Qa; to find 
that D''Q„6 = -|(aA:)-i(-A:2 + 3K)Qa, where as before D'^D^Q = -A^Q and A = | [8, 18]. 

^Notice that Do(lnTy ~ (Da lnT) +/f"(Da(lnr) + Aq) where Da = hab^^ and — !ia(lnT)- ~ Hiia- It is also 
well known that D''(Da InT) ~ (D^ InT) + HQ:)^ InT). 
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and taking another spatial covariant derivative, we obtain one of key equations in the primary 
source calculation: 

(D^ InT) + 2F(D2 InT) - ^B^D^'tc) ^ (D^Dbtr"'' - D„g°) - if(DM„) . (G.31) 



G.4. THE ADIABATIC CONDITION 

The entropy perturbation Sa, for a radiation-dust almost-Friedmann-Lemaitre universe, is 
given by 

iD,lnp^-iDalnp„ , (G.32) 

and this gives the relation used in [13]: 

D^^(,„~ (iD„lnp„-5„)p„ . (G.33) 

The adiabatic condition is then characterized by = ^. This can either be written as 
D„p„ = i?D„p„, ^STB-Qa = +i^D«Q„ and ST{k,t) = |A(fc,t). 

Provided the initial conditions produced, for example after inflation, lead to adiabatic pertur- 
bations, these perturbations will remain adiabatic until decoupling, however after decoupling 
generic density perturbations do not satisfy the adiabatic condition S(^rm) = 0. This is due to 
the fact that the average velocity of the radiation does not proceed along geodesies, while the 
matter does. Thus any perturbation that starts off adiabatic at last scattering will not remain 
so [13]. 



G.5. ALMOST-FRIEDMANN-LEMAiTRE SCALAR PERTURBATIONS 

In the case of scalar perturbations, the almost-Priedmann-Lemaitre Einstein field equations 
formatter domination reduce to the following CGI perturbation equations: 

&ab + 2Haab + Safe « <^ (aVafc)' « -O^Eab , (G.34) 
Kb + 2.HEab + \pM(^ab « , {pZ^Eab)' « —(^ab ■ (G.35) 

Taking the time derivative of the above equations, we obtain 

[a^f^ab) + Hia^Uab) « +^a^PM(^ab, (G.36) 

{p-'Eabr+2H{p-'Eab) « ^Kfe. (G.37) 

To see how these relate to evolution equations for the density gradient, we take the spatial 
divergence of (G.37) [8]. 

A useful consequence of the above relations is that for matter domination^, the monopole 
equations take on a simple form: 

(DalnT)- + H{Da\nT) + ^Dap^c) « -^oVafc . (G.38) 

*By adiabatic we mean that the comoving entropy density is constant. 

further consequence of (a^Baj)' « — ^a^pCj,;, and p = 3HqCIo<^~^ for 0,0 = +1; is that {a^E^i)' pa 
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To summarize, the matter dominated limit lead to the following simple relationships between 
the dynamical variables and the electric part of the Weyl tensor: 



Using 



CTab 
Da Inp 



0, 



(G.39) 
(G.40) 
(G.41) 
(G.42) 



E,, = J2^Qab = D{a^b)%{x) and D^x) =^^%{k,t)Qa (G.43) 



we obtain 



(G.44) 



Note that in [13], Eab = J2ki^^ /^^)^kQab which is based on the notation of Kodama and 
Sasaki [50]. The relationship between and the potential used here is $ = (A:^/a^)$fe. 
The evolution equation for the Newtonian like potential follows from (G.37) [26, 28]: 



(p-i$j + 2i/(p-i|.)«i|.. 



(G.45) 



In terms of the conformal time derivative {dt = adrj) which we denote by a prime '/' this 
becomes: 

(p-i$)" + i7(p-i$)'«ia#. (G.46) 
On rearranging (G.45) using (G.44) and (G.16-G.18) we find the evolution equation for $e: 



which can be simplified to give 



(a^yit, k) + 2H{a%) {t, k) « iH^VLoa-'%{t, k) 



(G.47) 



(G.48) 



Notice that for if = we obtain the usual equation $b ~ —AH^e for dust as used in [13]. 
This gives the well known result 



$,(/c,t)«*+(fc,0) + $-(fc,i)i-5/3. 

Considering the constant mode only, we then have : 

($p-i) « {3H^flo)-^a[k'^,{k,0)] « {3H^no)-'a[kHA{k,0)], 



(G.49) 



(G.50) 



It then follows that ($p~^)' ~ fc^d and ~ ~^'h ~ 0- Hence we recover the standard result 
that the 'potential fluctuations' $g are time independent for the K = matter dominated dust 
scenario. 

We can then write a^A{k,0) = D /a{a^A{k,0)) = {D / a)<^A{k,t) where D is the the linear 
growth factor. The flat dust case is recovered using D = a. 



{^p-^)^{3Hln,)-'^[kHA{k,t)\. 



(G.51) 
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APPENDIX H 
The Einstein Field Equations in the Energy Frame 

In the energy frame Qa = 0, one finds the following useful constraints and a simple form for 
the evolution equation for the electric part of the Weyl tensor: 

(diva)„ ~ +^D»e, (H.l) 

{divE)a ~ +^Ba{p + PR), (H.2) 
{pM+pR+p)Aa ^ -DaP, (H.3) 

These are very similar to those found in the case of matter domination [59] . The useful features 
that arises when using the energy frame are : Firstly, we have a simply relationship between 
the expansion perturbation and the shear (H.l). Secondly, the Newtonian like potential can 
be related to the density of the matter and radiation content (H.2). Thirdly, by taking the 
divergence of (H.3) we find that^ 

(pM + ^Pii)(divi)~-D% ^DVfl^-(l + i?-')p-'(divi) , (H.5) 

from which we can find an evolution equation for (div A) from the perturbation equations for 
the radiation energy density. Finally, equation (H.4) gives the relationship between the shear 
and the Newtonian potential. 

H.l. ON RELATING A^p TO THE CURVATURE AND Ct 

The Friedmann equation in a Friedmann-Lemaitre universe is 



H' + ^:^Ip ^ a^H' + K = ^a'p . (H.6) 
3 3 



Using definition of the density parameters: 

0=^,and^=^. (H.7) 

we can show that 

aV Sa'^H'^n . and a'^p, ~ 3a^H%. (H.8) 
Eq. (H.6) and eqn.(H.8) can be used to deduce that 

a^H^ + K = a^H^n,^a^H^ = — ^ . (H.9) 

[il 1) 



Using this and (H.8) we obtain 

(n-1) 



2 



^Note that this can also be written as D2(lnpfl) ~ |(1 + R){div A) using D^{lnpn) = 4D2(InT). 
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APPENDIX I 
The Correlation Function 

If {{SN/N)"^) is the square of the variance in the number of objects in a volume V, then the 
correlation function, the excess probability over a random variable of finding an object within 
a distance x of a given object [64], is given by 



^^^^ = w 



(1.1) 



where V{x) is the volume enclosed with a radius x ^or flat universe {K = 0). The power 
spectrum 'P(fc) is related ^ to the correlation function for a given distribution [49]: 

^(x) = ^/fc'*W^ ^ nk)=^^Jx'dxaxf^- (1.2) 



In an open universe the definition for the correlation function can still be retained. To see 
how, one first notices that in a flat universe {{SN/N)^)y(^^-^ {N)y^ ~ V~", where {N)y is 
the average number of objects in a volume V. In a space of constant negative curvature the 
volume enclosed by a sphere of radius x is Vix) = '''■(sinh(2x) — 2%) so that 



^(X) ~ V-'^ ~ (sinh(2x) - 2x)-". (1.3) 

The power spectrum for a power law correlation function in the volume in an open universe is 
then 

1 r> f . , 2 J sinfcx 1 



Vik) = —,B] -h^X^X^(-^^^J3^ , (1-4) 

where B is the normalization constant. This diverges for small x so a small scale cut-off 

is necessary [49]. To relate this to the power spectrum today on scales measured by galaxy 
surveys, the power spectrum is multiplied by the square of a transfer function T^{k). The 
power spectrum can then be normalized to o-g. 



APPENDIX J 

Power Spectrum Normalization 

There are two possible normalization schemes which one can follow. 

1. On considering a power spectrum of the form 

v{k) = Aikmr-' , (J.i) 

where tjq ~ 3to ~ '^Hq^ for f^o = 1 gives the conformal time today, the scale factor can be 
normalized and A is one way of expressing the amplitude of scalar perturbations since it is 
related to the dimensionless scale of matter fluctuations at horizon crossing, Xh, by [63, 2, 7, 80]: 

Xl = -A. (J.2) 

TT 

2. The alternative scheme is to consider a power spectrum of the form 



r{k) = Bk"T\k) , 

'^Ccire should be taken with the normahzation convention. 



(J.3) 



56 



GEBBIE, DUNSBY AND ELLIS 



where the transfer function T{k) ~ +1 on large scales. This means that if the fluctuations 
arise purely from the Sachs- Wolfe effect (potential fluctuations) near decoupling, B and A can 
be related at n = +1 [80] 

V{k) = 2tt^ rj^AkT^ik). (J.4) 

What is often used is V{k) = 2.5 x W^AT'^{k). The imits of A arc Mpc^ and of (Mpc^). 
For standard - CMD it is the convention to use the parametrized transfer function [6] : 

T{k) = 1 + (afc + {bk)i + {ckfY " , (J.5) 

where 

a = 6.4r-i/i-^Mpc , 6 = S.Or-i/i-^Mpc , c = IJT-^h-^Mpc , u = 1.13 . (J.6) 

Now the shape function F can be given as approximately F ~ Cloh. (by choosing h = 0.5 and 
flB = 0.05 the shape function is given as F = 0.48). 

Large scale flows provide a measure of the power spectrum in as much as the variance of the 
velocity field sphere of radius Xf, v^^g{xf) can be expressed as an integral over the power 
spectrum. On small scales (clusters of galaxies) the power spectrum is normalize to erg, the 
variance of the galaxy distribution on scales of Xf = 8/i~^Mpc [12]: 

^l = ^ = ^j k-'dkV{k)T\k)W\kxf) , (J.7) 

where bp is the Xf scale 'bias' such that cr™ = cr^/hp and the appropriate variance and 
W{x) = 3(sina; — x cos a;)l/a;'^ is the top-hat function. It should also be pointed out that it is 
more convenient to use the form 

.| = f f^(*™r«T^w(5^fe^). 

The variance of galaxies possibly biased to the matter {5gai = bdp) is roughly unity on the 
scale of 8h~^ Mpc. Interestingly enough the standard- CDM normalization from COBE seems 
to give CTg — 1.3 which is seems to imply that it is not correct to assume a pure Sachs- Wolfe-HZ 
power spectrum or even a n = 1.15 Sachs- Wolfe one. An appropriate table of standard-CDM 
normalizations is provided in [9]. 



APPENDIX K 
The Open Almost-Priedmann-Lemaitre Case 



K.l. EXTENDING THE INTEGRAL SOLUTION TO THE OPEN CASE 



In the main body of this paper we keep to the almost-flat-Friedmann-Lemaitrc case for 
clarity; the extension to the open case is straight forward. The point is that we shown that 
the standard results can be recovered from the CGI approach. Given that the generic linear 
Friedmann-Lemaitre models are well treated in the standard literature, we provide only an 
outline for the open case. The essence of the open solution is given via: 



4 + k 



1 



K 



{{£ + lf-l)]Te+i-Te-i 



{2£ + 3){2£+l) V 



(K.l) 
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The left hand side, the homogeneous case, is solved in very much that same way as for the 
flat case. 

The flat radial eigenfunctions are found from (K.3): 



(2£ + 1) 



(K.2) 



where Qim = 3e.{kx)yf^OAf For the general almost- Friedmann-Lemaitre (linearised 
Priedmann-Lemaitre) models [39, 35, 83] ^: 



d 



dx " (2^+1) 



(2£+l) 



1 - ^(^ + 2) 



(K.3) 



Thus, in order to construct the open solution, we just replace the homogeneous solution for 
z/^ + 1 = }p' I {—K) [35, 81] after reading off the solution using the recursion relation: 



(K.4) 



Of course one needs to be careful to redefine the wavcnumbcr. The integral solution (87) is 
then used in (K.l), i.e. , the solutions (K.4) are substituted into (93) to recover the open 
integral solutions. 

One carries out the same treatment as for if = but using (K.3) and (11-13). Alternatively 
one can re-deflne the mode expansion M^IQ] [18]. 

All that remains is to solve the evolution equations for the scalar perturbations, the coef- 
flcients Ci{r], k) now include curvature terms when written out in terms of the B terms (one 
uses the open recursion relation instead of the flat). In turn, these terms, B, will also pick up 
curvature terms when written out in terms of the perturbation variables. When the curvature 
starts to dominate the evolution, one gets an additional ISW contribution (which will be similar 
to the late ISW effect for a A dominated flat model). 

K.1.1. Extending the Power Spectrum to the Open Matter-dominated Ccise 

We consider the primary anisotropy term: 



(2^+1) - fc^i(^*)^aA7?*), 

Bi(r?B) ^ +^(afc)-i($(A;,77.)p~')(fc'-3if). 

This gives the Sachs- Wolfe effect in the open case. 
For spaces of constant negative curvature in an almost-Friedmann-Lemaitre model: 



(K.5) 
(K.6) 



($(r)$(fc'")) = (27r)^ 
Using (K.5-K.6) with if < we find 

^'^''°)(27TT)"" 
It follows that (noticing that Ajj* = x) 



K 



1^0 



2 3if - fc2 {^0 - 1) 



V{k). 



(K.7) 



(K.8) 



(2^+1)2 



(27r) 



fc2 (f^o-1). 



V{k). 



^Notice that (u^ + (i + + 1) = {1 - l{l -|- 2){K/k'^)). 



(K.9) 
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Then we can use 




(K.IO) 



for 



,2 



1-7:2 



(K.U) 



(2n+ l)(2n- 1) [ fc^ 



and (198) to find the angular correlation function. It is more straight-forward to redefine the 
mode function Gi[(5] in the mode function expansion, that is we use M^[(5] = {pi)^^ G i\Ci\ 
instead of G£[Q]. This then allows us to retain the flat like form for the angular correlation 
function (198), however we must retain the flat mode eigenfunction normalization: 



where = k^ — \. 

The angular power spectrum can then be normalized to a given structure formation theory 
such as standard-CDM. It would seem that the favoured model (by current observational limits) 
is the A CDM model [54]. 

Also, in the open case from (K.12 and 204): 




(K.12) 




(K.13) 



Here V* (206) is defined as before, while the mode expansion has been carried out in terms 
of rather than Ge[Q] [18]. 



